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ABSTRACT 


The flow field in which aircraft propellers must operate is spatially 
nonunifom due to propeller plane inclination and the presence of the 
airframe. Such a flow field alters the aerodynamic loading on the propel- 
ler blades and may cause excessive bladfe vibration, stresses, and acoustic 
noise. In order to design or analyze the performance of a propeller to 
be used on a given aircraft, the nonuniform flow field at the propeller 
plane must first be known. 

In this study, a computer program was developed to calculate the 
three-dimensional, steady, incompressible, inviscld, irrotational flow 
field at the propeller plane (propeller removed) located upstream of an 
arbitrary airframe geometry. The program uses a horseshoe vortex of 
known strength to model the wing. All other airframe surfaces are modeled 
by a network of source panels of unknown strength which is exposed to a 
uniform free stream and the wing- induced velocity field. By satsifying 
boundary conditions on each panel (.the Neumann problem) , relaxed boundary 
conditions being used on certain panels to simulate inlet inflow, the 

source strengths are determined. From the known source and wing vortex 

* 

strengths, the resulting velocity fields on the airframe surface and at 
the propeller plane are obtained. All program equations are derived in 
detail, and a brief description of the program structure is presented. 

A user's manual which fully documents the program is cited in the refer- 
ences. 

Computer predictions of the flow on the surface of a sphere and at 
a propeller plane upstream of the sphere are compared with the exact 
mathematical solutions. Agreement is good, and correct program operation 


is verified. 


iv 


Published experimental data are scarce. However, measured flow field 
data at a propeller plane of a twin-engine aircraft are present in the 
literature. Computer predictions for this aircraft are compared with the 
published test data. Reasonable agreement is observed, further validating 
the program. 

Results of a parametric study are presented which demonstrate wing- 
induced, aft fuselage-induced, and cowl inlet inf low- Induced effects on 
the flow field at the propeller plane of the single-engine Piper Cherokee 
PA-28-180 aircraft. Finally, a complete mapping of the computed flow 
field at the propeller plane of this aircraft is presented. However, no 
experimental data are available for comparison. 
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Chapter 1 / 

INTRODUCTION 

1.1 Statement of the Problem and Scope of the Investigation 

Improving the aerodynamic and vibrational characteristics of general 
aviation aircraft propellers has been an ongoing concern receiving re- 
newed emphasis in recent years. Nearby aircraft structures alter the 
airflow at the piano of the propeller. This nonuniform flow field in 
which the propeller operates produces fluctuating aerodynamic blade 
loadings which increase vibration and fatigue in the propeller structure. 
Also, the overall aerodynamic performance of the propeller may be reduced. 
Knowledge of the flow field at the plane of the propeller, as induced by 
the aircraft structure, enables the propeller designer to match the pro- 
peller with the airframe to reduce the severity of these difficulties. 
Thus, it is desirable to develop an analytic method for predicting the 
flow field in which the propeller will operate, for any arbitrary air- 
craft configuration. 

In reality, for a tractor propeller, there is a mutual interference 
between the propeller and airframe. The propeller slipstream produces 
a disturbed flow of higher velocity which impacts on the structures, 
nacelles or cowling, immediately behind the propeller. This situation 
produces extra pressure drag on the airframe which, considered alone, 
has a detrimental effect on the performance of the propeller-aircraft 
propulsion system. Simultaneously, however, the airframe components 
obstruct the flow entering the propeller disk to produce a nonuniform 
flow field at the propeller plane. The nonuniformity of this flow 
gives rise to increased propeller vibration. Flow blockage by the 
airframe can reduce the apparent advance ratio of the propeller and 
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improve the efficiency of the propeller. However, generally, the mutual 
interference tends to reduce the propulsive efficiency of the propeller- 
aircraft combination to an extent not immediately obvious. This mutual 
interference effect is very difficult to predict analytically and has 
usually, been investigated experimentally. 

However, as a starting point in exploring the propeller-aircraft 
interaction problem, this study confines itself only to examination of 
the influence of the airframe upon the flow field at the propeller with 
the propeller removed. It is this airframe interference which is dom- 
inant in creating the fluctuating aerodyanmic loads and vibration on the 
propeller blades. Propeller interference on the airframe and the mutual 
interference ace not considered. 

Only tractor propeller-airframe configurations are assumed in this 
study, as tractor propellers operating in front of the aircraft are more 
common than pusher propeller configurations. However, the methods pre- 
sented in this study are applicable to some pusher propeller configura- 
tions. 

A second assumption is that the flow is incompressible; thus, Maoh 
number effects are not introduced. 

Third, steady, inviscid, potential (irrotational) flow is assumed. 
This assumption is valid for propellers operating upstream of the airframe 
and not in the wake of any airframe components. Viscosity effects are 
insignificant in the flow at the forward regions of the airframe with 
the propeller removed. 

The flow field at the plane of the propeller is influenced by 
several factors. Two factors are aircraft angle of attack and the tilt 
of the propeller plane relative to the airframe. Many propellers are 
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mounted with some angle of attack and sideslip relative to the fuselage 
in order to improve stability. Angle of attack combined with propeller 
tilt causes the blades to see an azimuthal variation in velocity, and 
hence, fluctuating loads exist. This situation can exist even in the 
absence of any airframe interference. A third factor is the forward 
spacing of the propeller hub from the airframe as well as the vertical 
and lateral position of the hub relative to the airframe. The influence 
of the airframe on the propeller flow is lessened if the propeller is 
spaced farther from it. A last factor influencing the flow at the pro- 
peller plane is the shape of the individual airframe components and their 
positions relative to one another. 

Each airframe component has its own type of influence on the pro- 
peller plane flow field. The lifting wing upwash creates radial and 
azimuthal variations in the axial and tangential velocities at the pro- 
peller plane. Axial velocity tends to be increased over portions of the 
propeller plane located above the wing, while axial velocity is decreased 
over portions below the wing. Nonlifting components, the cowl-fuselage 
on single-engine configurations and nacelles on multi-engine configura- 
tions, tend to block the flow at the propeller. The axial velocity will 
be less than free stream velocity at the propeller plane, especially at 
the radial positions near the hub. Also, because these bodies are not 
usually bodies of revolution and because the propeller hub may not be 
centered in front of the cowling or nacelle, these bodies can also induce 
an uneven distribution of radial and tangential velocity components over 
the propeller plane. 

Combining the effects of all the above factors leads to a propeller 
plane flow field which is fully three-dimensional. The propeller plane 
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flow field cannot be simplified by assuming axisymraetry or assuming that 
aircraft components are bodies of revolution. 

The purpose of this study is to predict the three-dimensional, 
steady, incompressible potential flow at the plane of a propeller having 
any orientation and position in front of an arbitrary nacelle-fuselage- 
lifting wing combination. Development of a finite element numerical 
method is presented which calculates the flow on the surfaces of the 
fuselage and nacelles as well as at the plane of the propeller. This 
numerical approach is incorporated in a computer program which is des- 
cribed, in this thesis. 

Numerical predictions of surface pressure and velocity at the pro- 
peller plane are compared with exact analytical solutions i.vr a single 
body in order to check the method. 

Little usable experimental data is available in the literature for 
comparison. However, flow predictions at the propeller plane of a twin- 
engine aircraft are compared with the published wind tunnel test results 
for that aircraft. 

Computer predictions of the flow field at the propeller plane of a 
Piper Cherokee PA-28-180 aircraft are presented. However, no experi- 
mental flow field data for this aircraft are available for comparison with 
the predictions. 

A parametric study of the flow predictions for the Cherokee aircraft 
is made to determine the importance of the various airframe components 
on the flow field at the propeller plane. In particular, the acceptabi- 
lity of neglecting aft portions of the airframe distant from the pro- 
peller in making computer predictions is investigated. Also, the effect 



of wfng dihedral and inlet inflow at the cowl or nacelle upon the flow 
at the propeller is examined. 
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1.2 Previous Investigations 

Examination of the effects of afterbodies upon propellers has been 
done since the early 1920's. Most of the early work, prior to 1948, was 
both theoretical and experimental but was limited to body influences on 
the steady propeller thrust, torque, and propulsive efficiency rather 
than unsteady or vibratory loading. In much of this experimental work 
the actual flow field induced by the body at the propeller was not quan- 
titatively surveyed; rather, the propeller was operated in the flow field 
and its steady performance measured. 

Lesley and Woods (1) performed early wind tunnel tests of propeller- 
body interaction effects. Steady propeller thrust, torque, and propul- 
sive efficiency were measured with the propeller operating in front of 
disks and other flat-faced cylindrical afterbodies. Results indicated 
propulsive efficiency was less for the propeller operating with a slip- 
stream obstruction than for the free operating propeller. 

Durand (2) extended the work of Lesley and Woods Cl) hy use of 
actual airframe shapes in the propeller slipstream. A cowl-fuselage- 
high wing combination was tested at various spacings aft of the propel- 
ler. Increased propeller-airframe clearance reduced loss of propulsive 
efficiency, and the influence of the body fell rapidly as clearance 
increased. Durand concluded propeller-airframe interaction should bo. 
a design consideration. Also, use of fuselage shapes producing minimum 
interaction with the propeller is desirable, or if not possible, maximum 
propeller-airframe clearance should be used to reduce interaction. 
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Lock (3) first developed a two-dimensional analytic method for 
predicting propeller performance under the influence of a slipstream 
obstruction. The method was limited to a body shape consisting of a 
spheroidal nose attached to a semi-infinite cylinder aligned parallel 
to the free-stream velocity. Lock used the exact potential flow solu- 
tion for the body to get the radial distribution of axial velocity 
upstream in the propeller plane (propeller removed) . The axial velocity 
defect was incorporated in propeller blade strip theory to calculate 
the steady propeller performance with body influence. 

Lock (4) extended his own previous work. He analyzed the various 
power wastages by the propeller. To do this, a momentum analysis was 
used which incorporated the body- induced potential flow axial velocity 
at the propeller plane. Also, pressure drag on the body nose was used 
in the power wastage expression. This power loss expression could then 
be used with either propeller blade strip theory or vortex theory to 
predict steady performance with body influences. 

Lock (5) summarized his previous works. He applied his method to 
numerical examples involving spheroidal body shapes of various fineness 
ratios. He compared his predicted performance results with measured 
test results. 

Weick (6) made a limited survey of the axial velocity at the pro- 
peller plane of a Sperry Messenger aircraft (propeller removed) in a 
wind tunnel. At any particular test location in the propeller plane, 
the ratio of measured axial to free-stream velocity remained constant 
regardless of the wind tunnel test free-stream velocity. Weick concluded 
the variation of axial velocity radially f om hub to tip on the propeller 
plane to be significant and worth much consideration in propeller design. 
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Lesley and Reid (7) used a VE-7 aircraft fuselage and surveyed the 
distribution of axial velocity at the propeller plane (propeller removed). 
At each of three azimuth positions checked, the axial velocity varied 
from zero at the hub to near free stream at the propeller tip radius. 

An averaged radial distribution of axial velocity was presented. Next, 
several propellers of various pitch distribution were tested freely and 
then in front of the fuselage. Generally, each propeller operated less 
efficiently when in the presence of the fuselage. Significantly, one 
propeller having a blade pitch distribution shape similar to the shape 
of the measured radial distribution of axial velocity in the propeller 
plane suffered the smallest reduction of efficiency. Lesley and Reid 
concluded a propeller designed such that each blade section locally 
attains its optimum angle of attack, at conditions of maximum efficiency, 
is superior to the conventional constant speed propeller for operation 
in the presence of a slipstream obstruction. 

McHugh and Derring (8) tested several full-scale propellers opera- 
ting in front of a family of radial engine cowlings having various 
diameters. Effect of cowling to propeller diameter ratio on the steady 
propeller performance was investigated. Results were two-fold. First, 
for cowlings less than one-third the propeller diameter, propulsive 
efficiency did not vary significantly with variation in cowling diameter. 
However, as cowling size increased above one-third propeller diameter, 
propulsive efficiency rapidly decreased. Second, presence of a spinner 
on the cowling increased propulsive efficiency several percent above 
that obtained without the spinner. 

Stickle, Crigler, and Naiman (9) extensively tested three full-size 
propellers, each with Clark-Y blade sections, operating in front of 
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various cowlings. The five body shapes varied from a radial engine 
cowling without spinner to a very streamlined body having a large spin- 
ner. All bodies had circular cross sections and were at zero angle of 
attack. Flow surveys at the propeller plane of each body (propeller 
removed) yielded the radial distribution of axial velocity in which the 
propellers would operate. Powered tests of all propellers operating in 
these flow fields were made. 

For the bodies tested, it was concluded that the change in body 
drag due to the propeller slipstream negligibly influence propeller 
performance and could be disregarded. The remaining direct influence 
of the body-induced flow field on the propeller blade load' distribution 
was dominant. 

Their results led to three conclusions. First, the velocity field 
Induced by the body at the propeller plane must be considered when design- 
ing the propeller blade pitch distribution. Second, body^induced -flow' 
has the strongest influence on the inner blade sections. For round shank 
blades, reduced axial velocity due to the body alters the drag on the 
shank and affects efficiency. For blades with, shanks having airfoil shaped 
sections, the reduced velocity increases section angle of attack at the 
shank sections, so the body causes a change in blade load distribution. 
Third, efficiency gains obtained with. a spinner depend on -the 'velocity 
at the hub due to the body. The spinner is beneficial if hub velocity 
is high but not as helpful if hub velocity is low. 

Wing-induced effects on propeller loads were examined by Kuhn and 
Draper CIO) . A tapered wing semispan was tested with two propellers 
attached to negligibly small nacelles. When the outboard propeller op- 
erated at 90 to 100 percent chord ahead of the leading edge, the propeller 
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pitching moment was double that produced by the free propeller. However, 
no lateral propeller force was created. Thus, Kuhn and Draper concluded 
the wing upwash was not producing significant azimuthal variation in 
blade forces. Rather, the increased axial velocity induced by the wing 
over the upper half of the propeller plane and the reduced axial velocity 
over the lower half were found to cause the pitching moment increase. 

With the 1940's came studies examining the vibration of propellers 
due to oscillatory loads produced when operating in a nonunifarm .flaw * 
field induced by the wing. Such a study was conducted in Great Britain. 
In part I of the study, Pos.tlethwaite., Carter, Perring> and Diprose Q-H 
made theoretical predictions of vibrating modes in two-, three-, four-, 
and six-bladed propellers. These propellers were operating upstream in 
the calculated flow field induced by a wing having Joukowsky airfoil 
sections*. In part IT of the study, Forshaw.) Squire, and. Bigg Q.2]T tested 
a ':hree-bladed propeller operating in the nanuniform flow field created 
by blowing a narrow axial jet of air along one propeller plane azimuth 
location. Measured torsion-bending and whirling-bending vibration modes 
compared favorably with predictions. 

Later, Corson and Miller (13) examined the vibration of a pusher 
propeller whose blades periodically passed through, the wake of a forward 
body. The reduced axial velocity in the wake was introduced as an instan- 
taneous change of local blade section angle of attack in propeller blade 
element theory to predict unsteady loading. Actual blade forces were 
measured also. Blade vibratory stress was increased if the drag of the 
body (intensity of the wake) increased. Secondly, all else being equal, 
large diameter propellers experienced less wake-induced stress than 
smaller diameter propellers. 


10 


Most recently (1948 to 1953), the National Advisory Committee for 
Aeronautics published a series of studies dealing with detailed flow 
surveys at the plane of the propeller and these flow effects on the 
oscillatory loads and vibration in the propeller. 

Vogely (14) performed flight tests on a propeller inclined to the 
longitudinal axis of the fuselage. He measured radial and azimuthal 
blade load distribution and efficiency as affected by propeller thrust 
axis inclination. He compared theoretical predictions which did not 
account for any wing-airframe influences on the flow field, 

Vogely observed that the pure inclination of the propeller could 
produce oscillatory loads, and blade stress could be reduced, at a 
specific aircraft angle of attack, by attaching the propeller to the 
fuselage with the proper thrust axis inclination. Vogely pointed out 
that the airframe induces a significant spatially varying flow field 
in which the propeller operates which can incur further blade load alter- 
ations. Vogely concluded that detailed flow angle surveys at the plane 
of the propeller should be made to obtain accurate results in predicting 
propeller vibratory loads. 

Shortly thereafter, Roberts and Yaggy (15) performed the first very 
detailed survey of the flow field in the propeller plane (propeller 
removed) of a twin-engine aircraft. The full-size, unswept wing aircraft 
was tested in the Ames 40- by 80-foot wind tunnel. 

The test conditions run consisted of flaps extended, no flaps, and two 
different nacelle inlet velocity ratios. For each condition, a series 
of runs was made over a range of nacelle angles of attack from two to 
12 degrees. 
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Survey points at the left propeller plane were densely distributed 
over the full radius and at fifteen degree azimuth increments. Several 
flow quantities were measured at each point. Total local velocity was 
measured. The local velocity components were measured in terms of two 
angular quantities; angle of rotational flow was the ratio of local 
tangential to axLal velocity, and angle of outflow was the ratio of 
local radial to axial velocity. Also, at each point the upwash and side- 
wash angles, relative to the fuselage, were measured. Measured radial 
and azimuthal distributions of the flow quantities were presented. The 
distribution of upwash angle along the horizontal centerline of .the. pro- 
peller plane was plotted. 

Measured upwash along the horizontal centerline compared poorly 
with predicted upwash induced by the isolated wing modeled by lifting 
line theory. This comparison indicated the fuselage and nacelles con- 
tributed greatly to the upwash. 

Aerodynamic blade load predictions were made using steady blade 
element theory and included ’the. 'nonuniform flow field measurements in 
terms of Instantaneous spacially varying changes in local section angle 
of attack. Predicted loading varied approximately sinusoidally with 
azimuth, and when the measured flow field was used in the blade load 
predictions, those predicted load magnitudes were larger than loads 
predicted for an inclined propeller operating in free air. 

Nacelle inlet velocity ratio was found to influence the flow only 
at the central region of the propeller plane immediately upstream of 
the inlet. 

Roberts and Yaggy (15) concluded that the angle of rotational flow 
and, specifically, the upwash angle along the propeller plane horizontal 
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centerline were the most important quantities in producing the oscilla- 
tory blade loading. Also, they concluded it was desirable to develop a 
theoretical method for predicting the nacelle, fuselage, and wing con- 
tributions to the upwash at the propeller plane for any wing-fuselage- 
nacelle combination. 

Testing by Roberts and Yaggy (15) led to two further studies. Using 
the measured flow field data of Roberts and Yaggy (15) with steady pro- 
peller blade element theory, Rogallo, Roberts, and Oldaker (16) computed 
the airload variation for a propeller operating in the flow field. Based 
on airload predictions, the first order vibratory blade stresses were 
calculated. Blade stresses measured on the propeller operating on the 
twin-engine test aircraft compared favorably with predictions. This 
study gave two conclusions. First, steady state blade element theory 
was found adequate for predicting the magnitude and distribution of 
propeller oscillatory air loads provided the flow field induced by the 
body at the propeller plane was completely Known before hand. Second*., 
these predicted loadings allowed accurate prediction of first order 
vibratory blade stresses. 

Since Roberts and Yaggy (15) had verified that the upwash along 
the propeller plane horizontal centerline was the most dominant factor 
in inducing oscillatory blade loads, Yaggy (17) presented a theoretical 
method for predicting the upwash distribution. Yaggy’ s method was 
developed to predict the upwash contributions at the propeller plane 
due to the fuselage, nacelles, and wing of the test aircraft used in 
References 15 and 16. Lifting line theory was used for predicting 
wing-induced upwash. The nacelle was modeled as a simple closed body 
of revolution for which the exact potential flow could be calculated. 
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The portion of the fuselage in the vicinity of the propeller plane was 
modeled as an infinite elliptical cylinder for which an exact upwash 
expression could be derived from potential flow theory. Net predicted 
upwash at a point Ati the propeller plane was taken as the sum of the 
component contributions. Predicted upwash agreed fairly well with mea- 
sured results. 

Yaggy (17) indicated his method, derived for this specific aircraft, 
could be extended to other twin-engine aircraft of similar characteris- 
tics. The required characteristics were high wing aspect ratio, nacelles 
of circular cross section, fuselage cross sections resembling an ellipse 
at the propeller plane region, and that no two body components should 
be spaced closer than the larger body diameter. 

In three remaining publications, work was continued along the lines 
of the previous work (References 15, 16, and 17). This time the effort 
was to find the effects of wing sweep on the flow field induced at the 
propeller plane. Rogallo (18) and Rogallo and McCloud (19) extended the 
upwash prediction method of Yaggy (17) to the case of a highly swept 
wing-fuselage-nacelle combination. Rogallo and McCloud (19) also ob- 
tained limited experimental results for a swept wing aircraft which were 
compared with, predictions. Also, the swept wing results were compared 
with upwash measurements made with- a similar aircraft having unswept 
wings. 

Hoping to learn more about the effects of wing sweep on propeller 
vibratory loads, Rogallo and McCloud (20) performed detailed propeller 
plane flow surveys (propeller removed) with a swept wing-nacelle-fuselage 
combination. The test was as extensive as that done by Roberts and 
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Yaggy (15) . Wing sweep was found to induce first order vibratory propel- 
ler blade loads similar to those attained with an unswept wing. 

Since this work of the 1950's, there appears to be no published 
work presenting detailed propeller plane flow field surveys, or predic- 
tion methods. With the exception of Yaggy (17), Rogallo (18) , and Rogallo 
and McCloud (19) , whose prediction methods apply only to a small class 
of similar aircraft, there is no published method for predicting the 
nonuniform flow field at the propeller plane of an arbitrary wing- 
fuselage-nacelle combination. 


15 


Chapter 2 

THEORY OF THE POTENTIAL FLOW SOLUTION METHOD 
2,1 Historical Development 

With the development of high-speed digital computers in the 1950's 
and 1960's, it became possible to solve the potential flow about arbi- 
trarily shaped bodies using finite element techniques involving large 
systems of algebraic equations, 

Many such flow solution methods have been advanced since 1958 when 
Smith and Pierce (21) used source distributions on panels to model bodies 
of revolution. Hess (22) extended the method to bodies of revolution 
oriented perpendicular to the free stream flow, 

Hess and Smith (23) first solved the flow about nonlifting arbitrary 
three-dimensional bodies using distributions of sources over the body 
surface subdivided into discrete panel elements. Smith (24) explained 
the theory of Hess and Smith (23) and presented sample calculations. 

A thorough discussion of the panel element computer methods was done 
by Hess and Smith (25) , They discussed the theoretical details for pre- 
dicting both two-dimensional and fully three-dimensional potential flow. 

Throughout the 1960 's and 1970 ' s, many sophisticated computer pro- 
grams were written to predict the flow about three-dimensional lifting 
wing-body combinations. Several examples of such work are Hess (26), 
Woodward, Dvorak, and Geller (27), and Woodward (28). The theory be- 
hind these methods is the same, but each uses different singularities 
over the body panels. Sources, vortices, doublets, or combinations of 
these singularities have been used in attempts to improve the accuracy 
of the method. 
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All of such published computer program packages are limited to 
body surface velocity, pressure, force, and moment predictions. None 
of them have been written to calculate the flow at locations, such as 
at a propeller plane, away from the body surface. 

2 . 2 Theoretical Overview 

The method used here involves the solution of the Neumann problem 
for the potential flow around an arbitrary body-wing combination in a 
uniform free stream velocity field. Hess (25), Rubbert and Saaris (29), 
and Hess and Faulkner (30) explain the theory which is also summarized 
here. 

Steady, incompressible, inviscld flow is- assumed. Also, irrota- 
tlonality is assumed. Thus, the velocity field is the negative gradient 
of a scalar potential function. Three components comprise the velocity 
field. First is the onset free stream velocity, V, assumed to be. con- 
stant everywhere. Second is the velocity induced by the lifting wing; 

V v , if present. Third is the perturbation velocity induced by the body 
or bodies. 

A distribution of source strength density, a, over the body surface, 

vs 

S, corresponds to a scalar potential function, <)», which satisfies the 
Laplace equation and whose negative gradient is the velocity at any 
point. However, the wing surface is not included as part of the body 
surface. Rather, the wing is modeled by a specified vortex system re- 
sponsible for the wing- induced velocity. This velocity is thus a known 
quantity. 

At a point P on the body surface, the potential due to a unit 
strength point source at point q is 


l/[4irr(P,q)] 
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where r(P,q) is the distance from q to P. Then If surface S has a 
distribution of source strength density, o(q), the total potential at 
P, <f> (P) , due to the body or bodies is 


*(P) 


's 


q(q) 

4irr(P,q) 


dS 


( 1 ) 


where dS is an elemental body surface area. This potential also satis- 
fies the Laplace equation as a consequence of incompressibility: 


V 2 <KP) - 0 . (2) 

Also, due to the irrotationality condition, the perturbation velocity 
at P is given by -V<J> (P) . Thus, the continuity equation may be expressed 
as 


V • [~V<KP)] - 0 . (3) 

Equation (1) , which satisfies Equations (2) and (3) , is used with speci- 
fied boundary conditions on surface S to solve for the unknown o(q). 

The boundary conditions are that the normal velocity cn the body 
surface must be some prescribed value, F(q), while the potential, <j>, at 

A 

infinity must vanish. Generally, F(q) is zero at locations where S is 

a solid boum’.»ry but may be nonzero over areas of the body wall which 

are considered to be permeable. F(q) is the normal component at point 

q of the sum of the free stream velocity, V, the wing- induced velocity, 

V w (q) , and the perturbation velocity given by the negative gradient of 

Equation (1) . Taking the normal derivative of Equation (1) , using care 

in the integration at the singularity of P ■ q as explained by Hess (25) , 

and combining it with the normal component of V and V gives the follow- 

w 

ing boundary condition equation at point P on S: 

- i f + fy p » - F(p > ■ «> 


2 
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where n(P) Is the outward unit normal vector to surface § at P, and 3/3n 
denotes the derivative in the direction of n(P) . 

Solving the integral Equation (4) for the unknown o on the body 
surface concludes the solution of the complete flow field. 

Instead of analytically solving Equation (4) , which may be impossi- 
ble for complex surfaces S, it is solved numerically by dividing surface 
S into N discrete elements. Applying the surface boundary condition to 
each of the N elements approximates Equation (4) by a set of N simultan- 
eous linear algebraic equations for the N unknown values of cr. Solving 
for o on each surface element completes the flow solution, and the 
velocity anywhere on or away from the body surface may then be calculated. 
The details of this numerical method are explained in the remainder of 
this chapter. 


2 . 3 Body Surface Model and Free Stream Velocity Components 

All fuselage and nacelle surfaces are partitioned into a network of 
N plane triangular and quadrilateral panels. Figure 1 illustrates the 
panel network on a single body. 

S^ denotes the area of the ith panel. The panel control point is 
denoted by (XC^, YC^, ZC^) and outward unit normal vector by n^. Also, 
on the ith panel there is a constant distribution of source strength 
density, (a source strength-per-unit area) , see Figure 1. Appendix A 
provides the details of generating the panel network and obtaining panel 
areas, control points, and unit normal vectors. 

In terms of the body-fixed Cartesian coordinate system, the outward 
normal unit vector is written as 

n . ■ n ? + n 3" + n lc . 

i >v * t 


( 5 ) 
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As Is also explained In Appendix A, the body surface may be more 
simply modeled If all body cross sections are symmetric about a common 
plane of symmetry. In this situation the left (-y) side of the geometry 
is provided and the right side paneling is a mirror image of the left. 
Fewer panels need be considered, and this simplifies the problem, as 
will be explained in later sections of this chapter . 

On certain panels, say the ith one, inlet or outlet flow is modeled 
by specifying a value F^. F^ equals the amount of velocity allowed to 
pass normal to and through panel i. F^, as a positive quantity, repre- 
sents velocity inward through the panel. F^, as a negative quantity, 
represents outflow velocity. A solid boundary panel i has a zero value 
of F^ which is assumed and need not be specified. 

Referring to Figure 1, the body is Immersed in a uniform onset 
velocity, V, and is oriented in this velocity field with a specified 
angle of attack, a, and sideslip, g. 

Relative to the body-fixed axis system, the body, with orientation 
a and 8, senses three components of the onset free stream velocity ex- 
pressed vectorially as 

$ - ui + v j + wlc , (6) 

where the components are expressed as 

u ■ -V cos a cos ti , (7. a) 

v * +V cos a sin q , (7.b) 

and 


w * -V sin A , (7.c) 

where V is the free stream velocity magnitude, and Euler angle, n, is 
expressed as follows: 


q 


tan ^ (tan 2 cos o) 


( 8 ) 
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Appendix B derives the Equations (7) and (8) for the free stream 
velocity components. 

2.4 Body Panel Source- Induced Velocity Components 
2.4.1 General case of a non symmetric body panel network 

Expressions for velocity induced by the source distribution of the 
ith panel on Itself and by all other panels j on the ith panel are devel- 
oped as follows. Refer to Figure 1 and also Figure 2. On panel i, under 
consideration, the source strength density cr^ remains distributed over 
the panel in determining that panel’s self-induced velocity. However, 
in obtaining the velocity induced at panel i by all other panels, j, the 
source distribution is considered to be lumped into a point source at 
each of the remote panels j . 

As shown in Figure 2, the source strength distribution o^ produces 
only a velocity flux normal to panel i. This total velocity flux is 
numerically equal to the 

where V *ii is the normal 
the quantity of interest 

To compute velocity 
strength distribution on 

Qj , the strength of 

equals the volume flow rate of the source. On a sphere of radius r 


source strength density per unit area: 

a /S - 2V /S , 
i i n 1± i 


(9) 


velocity induced by panel i on itself. This is 
and is found to be from Equation (9) , 


ii 


ii 


( 10 ) 


induced at panel i by panel j , lump the source 
j into a point source of strength given by 


Q 3 ■ 


( 11 ) 


the point source in three-dimensional space. 
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with the source at the center of the sphere, the source-induced velocity 
is directed radially and is given using Equation (11) hy, 


v r “ (cr^/W 


( 12 ) 


equals the negative gradient of the potential, ^ , at radius r 
due to the point source. Conversely, ^ at radius r equals the negative 
of the line integral of from the source to radius r 


♦j * -/ V r dr . 


(13) 


Using Equation (12) in Equation (13) and integrating gives the 
expression for the potential at the control point of panel i due to the 
source at the control point of panel j as follows: 


a .-Zlfl 

*ij 47rr i;J ’ 


( 14 ) 


where r^ is the distance from control point j to control point i given 
by 

X ±i - £ (XC ± - XCj) 2 + (YC i - YCj) 2 + (ZC i - ZCj) 2 ] 1 ^ . (15) 

Then the velocity induced at panel i by panel j is the negative 
gradient of the potential as follows: 


-7* 

3x 3y 2 3 7 


(16) 


lj T lj 

Substitution of Equations (14) and (15) into Equation (16) gives 
O S 

v ±j - -J-ij[(XC i - XC.)t + (YC i - YCj)! + (Z^-ZCjtf] , (17) 


ij 


where r^ is given by Equation (15) . 

Equations (10) and (17) apply to the situation in which all N 
panels of the body must be considered individually which occurs when 
the panel network is nonsymmetric. 
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2.4.2 Special case of a symmetric body panel network 

In the special case of body symmetry, each panel on the left C-y) 
side of the body has a corresponding mirror image panel whose control 
point differs only by the sign of the y coordinate, see Figure 3. In 
this situation. Equations (10) and (17) are modified and applied only to 
the N/2 panels on the left side of the body. The modifications account 
for the induced effects of the right side "image” panels, and these 
modified expressions are uow derived. 

First, Equation (10), for panel self-induced velocity, is modified. 
In the symmetric body case, the velocity induced on left side panel i by 
itself, Equation (10), rust additionally include the normal component 
of velocity induced at panel i by a point source at the image of panel 


i ’ 1 image' 

Referring to Figure 3, the distance rjj. between panel i and 

image 

its image is 


ii 


- 2YC. 


image 


(18) 


Substitute Equation (18) into (17) and replace point (XC^, YC^ , 
ZCjj) by the control point on panel i^ ma g e to get the total velocity, 


i*- ima ge* induced by the point source at i 4 


on panel i. Then get 


image 

the component of velocity normal to panel i by taking the vector product 


of Y U < and n^using Equation (5) . The result is 


a i S i YC i n yi 


'ii 


image * 16 tt | YC^ j ‘ 


(19) 


Finally, to modify Equation (10) for the effects of the image 
panel, simply add Equation (19) to the original Equation (10) to give 
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°i S.YC.Hy. 

v ii " 2 [1 + 7I 3 ’ < 20 > 

in which panel i is on the left (-y) side of a symmetric panel network. 
Note, here is the velocity induced normal fio panel i by itself and 
its image panel. 

Now a modified form of Equation (17) for velocity induced by remote 
panels is derived for the case of body panel symmetry. Again refer to 
Figure 3. At real panel i there is a potential, $, due to the point source 
at real panel j, and there is also a potential at panel i due to the 
point source at image panel j. Due to symmetry, the source strength 
at panel j equals that at Also, note that YC^ ■ -TCj image’ Then 

denoting by r ij^ ma g e the distance from image panel j to real panel i, 
write: 


image 


[(X^ - XCj) 2 + (YC^ + YCj) 2 + (ZC t - ZCj) 2 ] 1/2 . (21) 


Applying Equation (14) to both panels j and j. , the total potential 

mage 

at panel i is. 


’ image 


ij 4irr. 


where r^ is given by Equation (15) , and tij.j^gg is given by Equation 


The velocity at panel i due to panels j and , is found 

by substituting Equations (15) , (21) , and (22) Into Equation (16) . This 


gives 
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0 4 S . (XC. - XC.) (XC. - XC .) 


L 


X 


ij 


(r 


ij 


)' 


rl i + I- 


(YC. - YC.) (YC. + YC.) 


A 


image 


r ij' 


ij 


image 


j 

r 


(ZC. - ZC.) (ZC - ZC ) 

+ [ i-—-L + i L-) k } . 


ij 


‘ij 


V 


(23) 


image 


This equation is the modified version of Equation (17) and is applied 
only to real panels i and j on the left side of a symmetric panel net-* 
work. 


2. 5 Wing Model Description and Wing- Induced Velocity 
2.5.1 Horseshoe vortex model and its induced velocity 

Should a lifting wing be present on the configuration , the surface 
of the wing is not modeled with source panels. Rather, it is represented 
simply by a single horseshoe vortex. 

Referring to Figure 4, the wing operates in the free stream velo- 
city, V. Wing span, b, root chord, c, and location of the quarter chord 
of the wing root, (X r , Y^, Z y ) , relative to the body-fixed coordinates 
are specified. Also, the quarter chord line may have some small angles 
of sweep, A, and dihedral, A. The horseshoe vortex representation has 
the bound vortex filaments attached to the quarter chord line. The two 
trailing vortices are spaced to give the horseshoe a span equal to irb/4. 
This spacing represents the mean span between the rolled up vortices 
trailing from an elliptic wing. These trailing vortices, which actually 
extend infinitely downstream, are truncated to a length, of 100 c for 
purposes of numerical modeling. 

Corresponding to the given body angle of attack, a, the wing lift 
coefficient, C. , is specified. C is assumed equal to the wing root 
section lift coefficient. Then the strength of the horseshoe vortex 
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filaments, r, is a constant which equals the total bound circulation at 
the wing root. Thus, using the Kutta-Joukowski theorem relating lift to 
bound vortex strength, the horseshoe vortex strength is given by 

cC.V 

r - -f- . ( 24 ) 

Note, r is a specified constant and not an unknown. This follows 
from the assumption used that the wing can influence the flow at the 
body panels, but the body will not influence the flow over the wing. 

At a point i of interest, the velocity induced by the wing horseshoe 
vortex is given in the body-fixed coordinate system as 

$ “ulf + v^ + wlt. (25) 

w i w i J 

The wing-induced velocity is calculated using the Biot-Savart law 
by applying the law to each of the individual straight vortex filaments, 
in turn, and summing the individual filament-induced velocities, W. Re- 
ferring to Figure 5, the Biot-Savart law for the velocity at point i loca- 
ted distance 'h- from a filament of length )t| is given vectorially as 

^ ’ <“> 

As shown in McCormick (31) , the magnitude of this velocity is given 
by 

\\\ * 4^< cos ®i + cos V ‘ (27) 

Should a vortex filament pass through or near the point of inter- 
est, distance h will be zero or very small causing numerical problems 
with the Biot-Savart law. Two measures are taken to prevent these prob- 
lems. First, if distance h is zero, the velocity components induced by 
the vortex filament are set equal to zero, and the Biot-Savart law. 
Equation (27) , is not used, thus preventing attempted division by zero. 
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Secondv an upper limit is imposed on the magnitude of velocity which 
a single vortex filament may induce. The upper velocity limit is 20 
percent of the free stream velocity. Thus, when distance h is such that 
Equation (27) yields a velocity magnitude greater than 0.2 V, the in- 
duced velocity magnitude |w^| is set equal to 0.2 V,. and velocity com- 
ponents are adjusted accordingly. 

Choice of 0.2 V for the velocity limit was made because this approx- 
imates the value of downwash velocity at the root three-quarter chord 
location of a wing with flat plate airfoil sections having an aspect 
ratio of 3 and operating at ■ 2.0, as calculated by thin airfoil 
theory. Thus, at body panels near the wing root region, where calcula- 
tions with Equation (27) will likely require the use of the velocity 
limit, 0.2 V is an appropriate velocity limit. 

Appendix C presents .expressions, employing the Blot-Savart law and 
the 0.2 V velocity limit rule used, by a computer subroutine to calcu- 
late the induced velocity at any point due to a straight vortex filament 
having any specified orientation. 

2.5.2 Justification for the use of the horseshoe vortex model 

To show the validity of using a simple horseshoe vortex model in- 
stead of more complicated vortex lattice or lifting line models, a thin 
airfoil will be investigated. 

A flat plate airfoil with the vortex attached at the quarter chord 
represents the two-dimensional case of a horseshoe vortex. Using Weis- 
singer’s approximation, the upwash velocity, jj, may be calculated. 
Analagous to a more exact finite wing model using a vortex lattice, the 
flat plate airfoil may be modeled using m segments, each with a vortex 
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at the segment quarter chord. Weissinger's approximation applied to each 

element can be used to solve- more accurately the upwash, y , ahead of the 

k in 

airfoil. 

Figure 6 compares the upwash at various upstream positions, xfc, of 
the single element airfoil to that for a multi-element airfoil having 
m ■ 100 segments. Comparing the approximate to the nearly exact model, 
there is a 13 percent difference in upwash at points at x/c ■ 0.4. Ty- 
pical twin engine aircraft propellers are positioned approximately at 
X/c ■ 1.0, where the upwash difference has dropped to 3.6 percent. Typi- 
cal single engine aircraft have propellers positioned approximately 1.5 
chords oir more ahead of the wing. As Figure 6 shows, the upwash differ- 
ence is 2.0 percent or less at those positions. 

As the main purpose of this study is to predict the flow at an up- 
stream propeller plane. Figure 6 Indicates that the use of the single 
bound vortex (horseshoe vortex) representation of the wing is justified. 


2.6 Boundary Condition Equations 

2.6.1 General case of a nonsymmetric body panel network 

Each point source potential, Equation (14) , on the body vanishes 
at infinity which satisfies one of two boundary conditions. 

The other boundary condition is applied to each panel on the body. 
This condition is that the normal component of velocity at the control 
point of each ith panel must vanish or equal a specified value of inflow 
velocity, F^, through the panel. The velocity at the panel is composed 
of free stream, wing-induced, and all N panel-induced velocity contri- 
butions. For the ith panel, the surface boundary condition is given as 


} . t + [ ? V • n 1 + V + V • n - {°* 80114 boundar y (28) 

w. n i ^j*l V ij i J ii i -Fji relaxed boundary’ ' K ' 

1 OH) 
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where is a positive value if it is inflow velocity and negative if 
outflow velocity. By substitution of Equations (5), (6), (10), (17), 
and (25) into Equation (28) , the boundary condition on panel i becomes 

+ £ jli ^< xc i - x Vv + < YC i - Y W< + « c i - z Vv> 

(Jfi) r ij 1 1 1 

■ -2[(u + u )n + (v + v )n + (w + w )n ] 
w i x i w i y i w i z i 

, ,0, solid boundary 

-2F,j» relaxed boundary " ■ ' 

By writing Equation (29) for each of the N body panels, a set of 
N simultaneous linear algebraic equations is formed in which the source 
strengths, a, are unknowns. In matrix notation^ this set is expressed as 

[A^HOj] - [B ;J ] i,j - 1, 2, 3, ... N . (30) 

Note, the [B] matrix containing the wing-induced velocity components 
and free-stream velocity components is the only quantity directly depen- 
dent upon the input values of a and 8 (Equations (7) and (8)) and C^. 

The matrix [A] consists of geometric parameters and is unaltered by a, 

8, and C T variations. 

A given input combination of a, 8, and will result in a particu- 
lar flow solution, [a] matrix. Assuming there are M total input combina- 
tions of a, 6, and C^, the system of boundary condition equations can be 
augmented to contain M sets of [B] coefficients and M sets of unknowns 
[ Cj 3 • Denoting by subscript k the kth input combination a^, 6^, and 
C^, Equation set (30) can be expressed In matrix notation ‘as 

tV'V ■ [ V • « l -*> 

where in accordance with Equation (29) the coefficients are 
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ij 


1, for (i - j); i, j - 1, 2, 3 ... N 


A ij" 


^tUXC. - XC.)n + (YC. - YC )n + (ZC. - ZC )n ] 

2^^ 1 3 x i 1 3 y i 1 3 *i 


for (i + j); i, j - 1, 2, 3 ... N , 


(31. b) 


and 


B ik ' - 2t( \ + + <v k + v »l k N 1 + ‘"k + "W 


+ ,0, solid boundary 

-2F^, relaxed boundary ’ 

for 1*1,2, 3 ... N ; and k * 1 , 2 , 3 . . . M 


(31. c) 


where is a positive value if it is inflow velocity and negative if 
outflow velocity with given by Equation (15); u^, v^, and w^ are 
direct functions of and 8^, by Equations (7) and (8). Also, uw^, v w ^» 
and w w ^ are the wing-induced components at panel 1 for the wing operating 
at a s per Section 2.5 of this chapter. 

Equation set. (31) represents the augmented set of N simultaneous 
equations to be solved for M different sets of N source strengths, a. 

With all a values known, the flow around the body will have been solved 
for each of the M input body orientations. 


2.6.2 Special case of a symmetric body panel network 

In the case of a symmetrically paneled body, only the N/2 panels on 
the left (-y) side of the configuration need be considered. In this case 
only a half-size system of N/2 simultaneous equations must be solved. 
Analagous to Equation set (31), this half-size system of boundary condi- 
tion equations is written by substituting Equations (5), (6), (20), (23), 
and (25) into Equation (28). However, summations are only made to N/2 
rather than to N. The half-size system in matrix form is 
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* A ij 3 *°jk 3 “ lB ik 3 * 


(32. a) 


where 


S i YC i n y 


ij 


[1 + — 4 ] , for (i - j); i - 1, 2, 3 ... N/2 , (32. b) 

o— I I 


Sit YC 


A ij - 2^[5^1 3 [« c t - + < IC 1 ' + ‘ ZO i ' 


+ ] 3 [(XC 1 - XC j )n x + (YC i + YC^n + (ZG i - ZC^n, ]> , 

-o , i i ^ 


^•^image' 

for (i + j); i, j - 1, 2, 3 ... N/2 , 


(32. c) 


and 
B 


ik ■ - 2[<u k + u "i k ) % 1 + (v k + V “l k )n ! r l + ( ” k + "W 


, rO, solid boundary - 

-2F^, relaxed boundary ’ k 


i - 1, 2, 3 ... N/2 and 
1, 2, 3...M , 


(32. d) 


where index i represents body panels on the left side of the configura- 
tion and where F^ is a positive value if it is inflow velocity and nega- 
tive if outflow velocity. Also, r^ is given by Equation (15) and 
r ij image 8*- ven by Equation (21). Components u^, v^, and w^ are func- 

tions of 01 ^ and 8^ by Equations (7) and (8). Also, u w ^, v w^» an< * w w^ 
are wing-induced components at left side panel i due to a wing opera- 
ting at as per Section 2.5 of this chapter. 

Equation set (32) is solved for M sets of unknown source strengths, 
o, on each of the N/2 left side body panels. The solutions on the right 
side image panels match those on the left side. The flow around the 
body will then have been solved for each of the M input body orienta- 


tions. 
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2.7 Body Surface Flow Field Velocity and Pressure 

Having solved the system of Equations (31) or (32) for the M sets 
of unknown pahel source strengths, a, the resulting flow velocity and 
pressure at each of the N panel control points is calculated using ex- 
pressions developed below. 

In the following, assume only one set of N solutions, o, has been 
solved. That is, set M equal to one in the equation sets (31) and (32). 

At the ith panel control point (XC^, YC^, ZC^), the resultant flow 
velocity vector, V b ^, is denoted by its body-fixed Cartesian axis com- 
ponents as 

V, - u. 1 + v. J + w. t . (33) 

b. b. b. b. 


This resultant velocity is the sum of the free stream, wing-induced, 
panel self-induced, and remote panel- induced velocities; 


V, - V + V + v.. + T V. . 
b. w. ii j-1 ij 

01 * 1 ) 


(34) 


The self-induced velocity, V^, is always directed normal to panel 
i. Thus, its Cartesian components are the scalar multiples of the unit 
normal vector at panel i. Using Equations (5) and (10) , is given 
in vector notation by 


'ii 


4 » t 


*i 


j + 


(35) 


Substitute Equations (6), (25), (35), and (17) for V, V w ^, V^, and 
respectively, in Equation (34) for the resultant velocity at 
panel i. Then according to Equation (33) the Cartesian components of 
V b± on the ith panel are 

V ■ + V> "TV* jJi -fe 1 [ H - 1 ■ (36 - a) 

UH) r ij 
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(v + v ) 

i 


N 

" j 3 1 

(j/i) 


s /j t YC i ~ tc j 

Air 1 3 


■ij 


(36. b) 


and 




■ Cw + w ) + 


w. 


n 


N 

Z 


dS) 


4ir 1 3 

ij 


(36. c) 


where u, v, and w are given by Equations (7) and (8) and where r^ is 
given by Equation (15) . 

The magnitude o£ the total surface velocity at panel i is 


'V ■ N + \ + \ 1 


2,1/2 


(37) 


Lastly, the pressure coefficient at the control point of panel i is 
given, using Equation (37) and free stream velocity, V, by 


C - 1.0 
Pi 


i j 2 


(38) 


Use of Equations (36) , (37) , and (38) on each of the N body panels 
completely defines the surface potential flow for a particular combina- 
tion of body. angle of attack, a, sideslip, 6, and wing C_ . 

Ij 


2.8 Flow Velocities and Flow Angles at the Propeller Plane 
2.8.1 Cartesian velocity components at a point on the propeller plane 
The propeller plane of radius R is centered at a point G^ub’ ^^ub’ 
Z hub^ relative to ttle body-fixed coordinate system. In general, the 
propeller plane is skewed with an angle of attack, o , and sideslip, & , 
relative to the body-fixed axes, see Figure 7. Although the propeller 
itself is not present, a right-hand sense of rotation in the propeller 
plane is assumed for velocity sign convention. 
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As indicated in Figure 7, a point P on the propeller plane has the 

specified radial and azimuthal position (r, i|») in a propeller plane 

cylindrical axis system. The cylindrical coordinates (r, t|/) of point P 

are transformed to corresponding body-fixed Cartesian axis coordinates 

(X , Y , Z ) using the following expressions which are fully derived in 
P P P 

Appendix D: 

X p m -r[sin tj* cos a p sin n p + cos <|) sin a p ] + \ab • <39 ' a) 

Y p - [r sin ^ cos n p ] + Y^ ub , (39. b) 

and 

Z p * r[sin sin a p sin n p - cos cos a p ] + , (39. c) 


where Euler angle, n p , is given by 

- 1 , 


n « tan [tan 3 cos a ] 
P P P 


(40) 


The total flow velocity vector, V , at point (X p , Y p , Z p ) is denoted 

by its body-fixed Cartesian components as 

$ * u t + v jf + w lc . (41) 

P P P P 


V is the vector sum of the free stream, wing- induced, and panel 
P 

source- induced velocities, V at the point P on the propeller plane; 

Pj 

V « V + V + ZV. . (42) 

P w p j-1 P3 

Equations (6) , (25), and (17) may be substituted for V, V , aad 

P 

V respectively, in Equation (42) if all i subscripts are replaced by 
PJ 

p and if point (XC ± , YC ± , ZC^ is replaced by (X p , Y , Z p ) . After making 
the aforementioned substitutions into Equation (42), the Cartesian com- 
ponents of the net flow velocity at point P on the propeller plane shown 
in Figure 7 are given as follows: 
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K S a X - XC 

u p ■ u + u w + jh «n [ 3 3 1 

p r pj 


(43. a) 


v - v 


m S.o, Y - YC. 

+ v„ + .£ -fc k P - J 1 


w p j-1 


pj 


(43. b) 


and 


N S 4°4 Z D “ 

w - w + w + Z -f-J-t- 2 =-J-] 

w p j-1 ^ T 3 

Pj 


P 


(43. c) 


where u, v, and w are given by Equations (7) and (8), and 

'pj “ [CX p " XC j )2 + (Y P " YC j )2 + (Z p " ZC j ) 2 l 1/2 • < 43 * d > 


2.8.2 Axial, radial, and tangential velocities at a point on the 
propeller plane 

At point P on the propeller plane, the Cartesian velocity compon- 
ents 0 ^, v p , and w p of Equations (43) may, by the following transforma- 
tions, be converted into axial, radial, and tangential components, v ap , 
v Ip , and v tp , respectively: 

v ■ [cos a cos n Ju + [sin n ]v » [sin a cos n ]w , (44. a) 

a p PPP PP PPP 

v » -[sin o cos ii> + cos a sin n sini|»]u + [cos n sin^lv. 
r p p P P P P P 

+ [sin a sin r\ sin - cos a cos i|»]w , (44. b) 

P P P P 

and 


v * [sin ot sin i|/ - cos a sin n cos i|>]u + [cos n cos 

t p p P P P P P 

+ [cos a sin \p + sin a sin n cos ^]w , 

P PPP 


(44. c) 


where n is given by Equation (40) , and if> is the azimuth, position angle 
P 

of point P. 


As indicated in Figure 7, v 3p is directed positive in the thrust 
direction perpendicular to the propeller plane. v tp is directed positive 
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radially outward in the propeller plane. v tp is directed in-plane per- 
pendicular to v r p and positive when directed clockwise around the pro- 
peller plane hub. Appendix D gives the complete derivation of trans- 
formation Equations (44) . 


2.8.3 Local flow angles at a point on the propeller plane 

The velocity components at point P on the propeller plane may be 

expressed in terms of four flow angles. These are angles of upwash, 

sidewash, outflow, and rotational flow denoted by e, t, X, and 6, 

respectively, and are functions of u , v , and w given by 

P P P 



P 


and 

i ” v 

t ■ tan” - ^-^-^] (46) 

U P 

such that e is positive when the flow is directed locally upward toward 
the negative z axis, and x is positive when the flow is directly locally 
to the right toward the positive y axis. 

The angles of outflow and rotational flow at point P are measured 
with respect to the propeller plane and were first defined by Roberts 
and Yaggy (15) . These angles are functions of v a , v r , and v t given, 

r r r 

respectively, by 


X - tan -1 [—^4 
a 

P 


^ W 

6 » tan 

a 

P 


(.47) 


C48) 


and 
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As shown in Figure 7, angle of outflow, X, lies in a plane contain- 
ing the thrust axis and radius line to the point P. X is positive when 
the lochl flow has a component directed radially outward. Angle of ro- 
tational flow, 6, Indicates the apparent tangential or rotational sense 
of the local flow at the propeller plane. 6 is measured from a line 
parallel to the thrust axis and lies in the plane perpendicular to the 
radius line to the point. That is, 6 is an angle which would appear in 
the propeller blade section diagram at that point. For sign convention, 
a right-hand propeller rotation is assumed, so a positive value of 0 
corresponds to an effective decrease in local blade section angle of 
attack from that encountered by the propeller operating freely in a 
uniform flow. A negative value of 0 corresponds to an increase in the 
local section angle of attack. Figure 8 illustrates the effect of 0 on 
the angle of attack of a blade section at radius, r, on a propeller rota- 
ting at an angular velocity of to. 

2,9 Three-Dimensional Flow Computer Program Description 

A computer program has been written which solves the subsonic po- 
tential flow about any arbitrary three-dimensional body-wing combination 
using the theory presented in this chapter. The surface geometry of the 
body, or group of bodies, must be provided in an organized set of dis- 
crete points. Geometry of the wing, if present, as well as the location 
and orientation of a propeller must be specified. Then for each, speci- 
fied combination of body angle of attack, sideslip, and wing C^, the 
program calculates the velocity and pressure coefficient at points on 
the body surface. Lastly, the program calculates the flow velocities 
and flow angles at each point in a mesh of points on the plane of the 
propeller . 
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Figure 9 presents the organizational structure of the program. It 
consists of fourteen subroutines each handling a major portion of the 
problem and linked by a main calling program. The program is coded in 
Fortran for use on the IBM 370/OS 360 computer system at The Pennsylvania 
State University. 

The program has been dimensioned to various sizes, the largest of 
which is capable of handling N ■ 2596 body panels and M « 6 combinations 
of a, 6, and C^, 

Several features have been incorporated. First is the symmetric 
input option feature (NSYMET ■ 0) which requires the input of only the 
left side of a symmetric body geometry. Image geometry is computed 
automatically and the program must solve a problem only half as large as 
would occur with a nonsymmetric configuration. 

Second, due to the massive amounts of input surface geometry 
required, a geometry check run feature (NCALC ■ 1) has been included. 

This feature is used when a new body geometry is being tried. The pro- 
gram generates and prints the panel geometry network which, allows the 
user to check for errors in paneling data and allows the user to identify 
the sequence numbers of panels to be characterized as inflow and outflow 
panels. Flow calculations, which, consume most of the computer time, are 
not made during a geometry check run thus reducing the amount of com- 
puting time which would have been wasted if flow calculations had been 
attempted with geometry data containing errors. 

Third r a feature (NPUNCH ■ 1) allows punched output to he produced. 
This feature, if selected, produces a punched card deck containing 
identifying information and propeller plane flow field data including 
axial and tangential velocity components at each, point on the propeller 
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plane. These output cards are formatted for direct application as Input 
data cards In the propeller performance analysis program of Aljabri (32). 

Because of the vast amounts of storage required by the program, It 
Is necessary to store the system of boundary condition equations, Equa- 
tions (31) or (32), on an auxiliary sequential scratch disk file. A 
second auxiliary file is also needed for printing certain output data. 

As a result, the system of equations cannot be solved by matrix inver- 
sion or other direct methods, but instead requires the use of iterative 
methods. The method of Gauss-Siedel iteration is used in this program. 

Additionally, because of the use of auxiliary files, slower program 
execution occurs. Thus,J.t'has been found useful to compile the program 
into a. highly efficient machine language deck using the Fortran H, optimi- 
zation level ■ 2, compiler available at The Pennsylvania State University. 
Actual jobs are performed by running this compiled machine language pro- 
gram. Execution time conserved by this technique can often mean the 
difference between success and failure in completion of the flow solution. 

The purpose of each, subroutine is given very briefly helow. 

INPUT reads and stores all input data. It also checks for certain 
input errors. 

VCOMP calculates the set of free stream velocity components using 
Equations (7) and (8) for each input set of a and g. 

EULER calculates angles- n, and rip using Equations (8) and (401, 
respectively. 

WGEOM generates the wing horseshoe vortex geometry based on input 
wing geometry if a wing is present in the configuration, 

PANEL generates the body panel network. It calculates all panel 
areas, control points, and unit normal vectors using the methods and 
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formulas given in Appendix A. It assigns a sequence number to each 
panel and counts the total number of panels. PANEL also generates image 
panels and establishes a system for matching "real" and "image" panels 
on symmetric configurations. 

COFSYM generates the matrices of coefficients, [A] and [B], for the 
half-size system of boundary condition equations, Equation set (32), 
used in the case of body symmetry. 

SOLSYM solves Equation set (32) for M sets of N/2 unknown panel 
source strengths using Gauss-Sledel iteration. 

COEFIC generates the matrices of coefficients, IA] and [B], for the 
full-size system of boundary condition equations. Equation set (31), 
used when the configuration is not symmetrically paneled. 

SOLVE solves Equation set (31) for M sets of N unknown panel source 
strengths using Gauss-Siedel iteration. 

WINGV calculates the velocity components induced at a point by- the 
wing horseshoe vortex. It uses subroutine VORTEX, below. 

VORTEX calculates the magnitude and components of velocity induced 
at a point by a straight vortex filament. Formulas given in Appendix C 
are used. 

VELOCI calculates and prints the surface flow velocity and pressure 
coefficient at each, panel control, point using Equations (361, (371, and 
(38). 

VPROPS calculates the propeller plane flow field. It first calcu- 
lates and prints u p , v p , and w p using Equation set (£3) at each, point on 
the propeller plane. Second’, it calculates and prints v ap , v Tp , and 

v tp using Equation set (.44) at each, point. Also, it calculates, using 

* 

subroutine ANGLES, and prints- flow angles c, t, X, and 6 using Equations 
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(A5) through (AS) at each point, Third, if desired, it punches values 
v ap» v tp> a **d other data on output cards for each point, 

ANGLES calculates flow angles e, t, X, and 0 for each point on the 
propeller plane using Equations (A5) through (A8) . 

A user's manual, Jumper (33), has been written which gives a com- 
plete description of program input and output. Also, the manual presents 
all necessary operating instructions and contains a complete program 
listing. 

However, the user's manual contains an older version of subroutine 
VORTEX which incorporates a quantity called h ||<n and which does not use 
the computational formulas given in Appendix C. 

Subroutine VORTEX listed in Jumper (33) should be modified by 
deleting the expressions and logic pertaining to h^ and then reorgan- 
izing the expressions to conform with the equations given in Appendix C, 

An errata sheet to the user's manual has been written which, describes 
all changes necessary to make expressions in subroutine VORTEX, originally 
printed in the manual, conform with the formulas found in Appendix C of 
this thesis. 
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Chapter 3 

TEST CASES OF THE FLOW ABOUT A SPHERE AND THE FLOW 
ABOUT A SIMPLE FUSELAGE SHAPE 


3.1 Exact Solution for the Potential Flow About a Sphere 

To ensure the three-dimensional potential flow computer program 
functions properly, the test case of the flow about a sphere has been 
used for comparison, because the exact solution for this flow exists. 

Figure 10 shows the cross sectional view of a sphere of radius, a, 
in a uniform velocity field, V. In terms of the polar coordinates, the 
velocity potential for this flow may be derived by mathematical solution 
of the Laplace equation. Also, the potential is found in Milne-Thomson 
(34) and is given by 

__ 3 

<|> “ V[r + — ■=■] cos y . (49.) 

8 2r 

On the sphere surface, only the tangential velocity component exists and 
is found, using Equation (49.), to be 


9<j> 


v 

Y 


If J « 7-a - I V sl " * ' 


(50) 


Then, using Equation (50) in the definition of pressure coefficient. 
Equation (38) , the sphere surface pressure coefficient distribution is 

C p - 1 - |lsin y] 2 • (51) 


The velocity potential may be rewritten in terms of the body-fixed 
Cartesian axes shown in Figure IQ. The potential becomes 

3 


4>s 


V[1 + 


2[ (x + a) 


2 + z 2j3/2^ * a * * 


(52) 
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To get the Cartesian velocity components anywhere on the x-z plane in 
Figure 10, the negative gradient of Equation (52) is taken. These x, y, 
and z components of velocity are given, respectively, by 


za lx + aj 




2[<x + a) 2 + z 2 ] 


2,5/2 


- 1]V , 


v - 0 
s 


(53. a) 
(53. b) 


and 


3 r 

v . "i 1 


a (x -f a)z 

r/ \2 . 2,5/2 

[(x + a) + z ] 


]V 


(53. c) 


3.2 Sphere Surface Pressure Prediction Compared With Exact Solution 
A sphere was modeled using ten equally spaced cross sections to 

generate 120 panels, and the flow was solved using the three-dimensional 
flow program. Figure 11 compares the predicted and exact, Equation (51), 
pressure distributions. Agreement is excellent over the entire surface 
with only a minute deviation near y * 90 degrees. This result verifies 
that the computer program functions properly. 

3.3 Flow Predictions at a Propeller Plane Upstream of a Sphere 

In addition to surface flow predictions, the flow velocities and 
Jingles were predicted at untilted propeller planes centered in front of 
the sphere. Results were predicted for two propeller plane-sphere 
spacings. Figure 12 gives the axial velocity distribution along the 
upper centerline of the propeller plane for both spacings tested and 
compares the results with the exact solutions. Equation (53. a). Figure 
13 presents the distribution of vertical velocity compared with the 
exact solution of Equation (53. c). Additionally, this figure combines 
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the axial and vertical velocity distributions by shtmng the corresponding 
upwash angle distributions. 

The results in both Figures 12 and 13 show excellent agreement between 
predictions and exact theory for both spacings at radial positions beyond 
0.6 sphere radii. At positions inside this radius, the predictions de- 
viate from the exact results. The deviation in upwash angle is the most 
pronounced. It is seen that as the propeller plane-sphere spacing in- 
creases, the computer predictions improve along the entire radius. 

At the inner regions of the propeller plane, the sphere surface is 
comparatively close to the plane, and the flow is dominated by the influ- 
ence of the nearby front ring of body panels. Jo. points at the inner 
radii of the propeller plane, which are near the body surface, the body 
panel network appears rather coarse and does not present as accurate a 
representation of the true body surface shape as is presented to points 
more distant from the body. This results in deviations of the predicted 
flow from the exact solution , at_ the inner radii of the propeller plane. 
However, at points at larger radii on the propeller plane, local panel 
domination diminishes, and the overall panel network more closely resem- 
bles a sphere. Thus, the computer predictions improve at points- farther 
from the body surface. 

Thus, it appears that in using the computer program for propeller 
plane flow predictions, the flow is most sensitive to the body surfaces' 
nearest to the propeller plane, and care must be taken to panel these 
surface regions densely and accurately. 
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3 . 4 Flow Predictions Upstream of a Simple Fuselage Shape and Effects of 

Ignoring Remote Body Surfaces 

Because the three-dimensional flow computer program requires large 
amounts of tedious body panel geometry input, it is desirable to reduce 
as much as possible the amount of this input without sacrificing accuracy 
in the propeller plane flow predictions. 

To examine the effect of neglecting remote regions of the body in 
predicting the flow at the propeller plane, a simplistic fuselage shape 
was used. This fuselage profile is shown in exact scale proportions In 
Figure 14. The cowl and tailcone have circular cross sections, while 
the cabin region has "pear-shaped" cross sections. A propeller plane 
is untilted and centered in front of the cowl with a spacing of 7,8 per- 
cent of the cowl length, f. 

Figure 14 compares axial velocity distributions at the propeller 
plane predicted by the computer program with the complete fuselage mod- 
eled and with cabin and tail removed. As might be expected, with the 
cabin and tail removed there is less flow obstruction, and a slight 
increase in the magnitude of axial velocity occurs. The velocity; distri- 
bution with the isolated cowling varies as much as 1.6 percent from the 
distribution predicted with the complete fuselage. This change in axial 
velocity considered alone might seem unacceptable. However, as the 
axial velocity changes so does the radial velocity component. Thus, as 
Figure 15 shows, the flow angularity, presented in terms of upwash angle, 
at the propeller plane changes only slightly by removal of the aft por- 
tions of the fuselage. Flow angularity shows the combined effect of all 
velocity components and can be considered a more useful parameter for 
describing the flow field as it would influence a propeller, 
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Results in Figure 15 thus indicate that in modeling a fuselage for 
propeller plane flow predictions, the body regions remote from the pro- 
peller, such as the tail and cabin, can be ignored without seriously 
affecting flow predictions. 

Paneling only the cowling of a typical single-engine aircraft fuse- 
lage should be sufficiently accurate for making flow predictions at the 
propeller plane. This will greatly simplify computer program input and 
decrease run time and cost. 

The effect of neglecting aft portions of a fuselage of an actual 
aircraft will be shown later in Chapter 5. 
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Chapter 4 

MEASURED AND PREDICTED FLOW FIELD RESULTS AT THE LEFT PROPELLER 
PLANE OF A TWIN-ENGINE AIRCRAFT 

4.1 Aircraft Geometry, Body Paneling, and Operating Conditions 

Roberts and Yaggy (15) conducted extensive wind tunnel tests in which 
the flow velocities and flow angles were measured at the left propeller 
plane (propeller removed) of the twin-engine aircraft shown in Figure 16. 
Because the experimental data was presented in detail, flow predictions 
for this aircraft were made using the three-dimensional potential flow 
computer program. Comparisons of computer predictions with the published 
experimental results were conducted to further check the accuracy of the 
computer program. 

Figure 16 shows a three-view drawing of the aircraft including over- 
all dimensions and the orientation of the left propeller plane. Figures 
17a and 17b present in detail the cross section geometry of the fuselage 
and nacelle surfaces located foreward of the wing leading edge. 

To generate the body paneling input data required by the computer 
program, the data in Figures 17a and 17b was used to obtain the paneling 
accuracy needed at the foreward body surfaces, which are in proximity to 
the propeller plane. At regions aft of the wing leading edge, which are 
more distant from the propeller plane, the surface geometry at a given 
longitudinal position was estimated by taking the cross section shapes 
shown in Figures 17a and 17b and scaling them to the dimensions indicated 
by the views in Figure 16. This produced a body paneling network of 
sufficient detail in the aft regions. 

Initially, the computer program was run using the wing vortex model 
and body paneling over the entire surface of the fuselage and nacelles. 
Unfortunately, with the body surfaces completely paneled, the program 
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failed to converge. Apparently, the matrix of boundary condition equa- 
tions generated by the complete Body panel network was not diagonally 
dominant or for some other reason was not solvable by the Gauss-Siedel 
iterative method employed by the program. Thus, it was necessary to 
reduce the complexity of the paneling network by excluding certain .por- 
tions of the aircraft geometry. 

Finally, after several trails, a simplified body panel network was 
devised which gave successful computer solution convergence. The wing 
vortex geometry was modeled as before. However, the right nacelle and 
aft portion of the fuselage were eliminated from the paneling network. 

Only the left nacelle and the fuselage nose, the crosshatched region in 
Figure 16, remained paneled. As previous results in Chapter 3 indicate, 
ignoring the right nacelle and aft fuselage, which are remote from the 
left propeller, should not significantly diminish the accuracy of the 
flow field predictions at the left propeller plane. The simplified left 
nacelle-fuselage nose paneling network consisted of 424 panels and re- 
quired the use of the nonsymmetric input option of the program. 

To match the aircraft flight conditions of Roberts and Yaggy (15), 
the computer model required a simulation of the engine cooling air flow 
through the left nacelle. All panels covering the nacelle inlet face were 
assigned a known constant inflow velocity, F, normal to and passing 
through the panels thus simulating inflow of engine cooling air. However, 
no cooling air outflow panels were specified. To obey the law of con- 
tinuity, it was assumed that the cooling air was exhausted infinitely far 
downstream or at the extreme aft location on the nacelle such that the 
exhaust air influence on the flow field at the propeller plane would be 
negligible. Based on this assumption, the presence of exhaust air could 
be ignored; thus, it was not simulated. 
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Table 1 summarizes the various flight conditions modeled using the 
computer program. Flight conditions 1, 2, and 3 on the table were also 
tested by Roberts and Yaggy (15). Flight conditions 4, 5, and 6 used a 
value of nacelle inlet inflow velocity which was not wind tunnel tested, 
and these conditions were run on the computer to investigate the effects 
of the inflow velocity parameter on the flow at the propeller plane. 

Table 1 also lists the geometric characteristics which remained con- 
stant for all flight conditions, including sideslip angle 8 , propeller plane 
reference radius R, and propeller plane orientation angles, a p and 8 p . 

The installed propeller sideslip angle, 8^, is zero for this air- 
craft. However, there existed a confusion about the value of the in- 
stalled angle of attack, a^, of the propeller plane. Roberts and Yaggy 
(15) conducted their experiments in the belief the propeller plane was 
not tilted downward with respect to the fuselage; that is. oip was believed 
to be zero. All of the experimental data was presented as a function of 
Ug, the geometric angle of attack of the propeller plane from the free 
stream, which under the assumption of zero a would also equal the air- 
craft angle of attack. However, Yaggy (17) later reported that following 
test publications' by Roberts and himself’ (15) it was found the prppeller plane 
had actually been oriented with a downward tilt of two degrees (a ■ -2 
degrees) during the testing. Discovery of this fact did not alter the 
measured results but produced an error in the presentation of the data 
in Roberts and Yaggy (15). In the published test results, data presented 
for a specified propeller plane angle of attack of a^, believed to equal 
the aircraft angle of attack a, are actually data corresponding to a pro- 
peller plane angle of attack of a_ but at an aircraft angle of attack 
equal to + 2 degrees. Similarly the published values tested cor- 


49 


respond to the aircraft angle of attack of o G + 2 degrees. Therefore, 
as Table 1 indicates, to model a wind tunnel test performed at a published 
propeller plane angle of attack of Cg, it was necessary to perform the 
corresponding computer run using an input aircraft angle of attack of 
a * Og + 2.0 degrees. 

4.2 Comparisons of Computer-Predicted and Experimentally Measured Flow 

Quantities at the Left Propeller Plane 

Figure 18 presents azimuthal distributions of axial velocity at the 
75 percent radius position for two extremes of body angle of attack. For 
the angle of attack of two degrees, there is fair agreement between 
computations and experiment. Agreement is much better for the 12 degree 
angle of attack case particularly at azimuths near the vertical centerline 
of the propeller plane. For both angles of attack, the computed distri- 
bution shapes agree well with the measured shapes. 

Figure 19 presents azimuthal distributions of flow angularity at the 
75 percent radius position corresponding to the velocities in Figure 18. 
Additionally, Figure 19 contains a set of curves for a midrange angle of 
attack of six degrees. Again the agreement between calculated and mea- 
sured distributions is good and improves with increasing angle of attack. 
The ability of the program to accurately predict the flow for higher 
angles of attack is good in that propeller vibration, which is more severe 
at higher angles of attack, may be more accurately calculated using the 
predicted flow field. Generally, in Figure 19, the flow angularity pre- 
dictions differ by approximately only one degree from the measured values, 
and at most only a three degree difference is seen for an angle of attack 
of two degrees. 
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As explained by Roberts and Yaggy 0-5) , the extremes in flow angu- 
larity occur near the horizontal centerline (near the 90 and 270 degree 
azimuths) of the propeller plane. The flow field near these azimuths 
produces the peak vibratory loading on the propeller blades. Additionally , 
along the inboard horizontal centerline of the left propeller plane, in- 
creased flow angularity is encountered due to the increased wing upwash 
in that region as well as due to the influence of the fuselage ‘nose near 
the inboard tip of the propeller plane. Thus, the ability to predict the 
flow along the horizontal centerline is of value for anticipating vibra- 
tion problems of propellers operating in the flow field. 

Radial distributions of angle of rotational flow along the inboard 
horizontal centerline (90 degree azimuth) of the left propeller plane 
are shown in Figure 20. Over the outer radii. beyond 0.25R, where a 
propeller blade would be more heavily loaded, computed flow angularity 
distributions agree well with the measured distributions particularly at 
the lower angles of attack. Increased flow angularity induced near the 
tip radius by the fuselage-nose may be seen on this figure and has been 
accurately predicted for the angle of attack of 12 degrees. However, 
computed and measured distributions diverge at radii less than the nacelle 
inlet radius in Figure 20. Fortunately, any weakness in prediction 
ability at these inner radii is not extremely bothersome, because a 
propeller blade is lightly or negligibly loaded at these radii. Propeller 
performance or loading calculations for blades operating in the predicted 
flow field of Figure 20 should not be seriously affected by some inac- 
curacy in flow field predictions at the inner radii. Disagreement between 
measured and computed flow angularity in the propeller plane hub region 
may be due to some shortcomings in the inlet inflow velocity modeling 
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technique used by the computer program. Inlet Inflow velocity effects 
are examined later in this chapter. 

Figure 21 shows the axial velocity distributions along the Inboard 
horizontal centerline for angles of attack of two and 12 degrees corres- 
ponding to the flow angles of Figure 20. The predicted and measured 
velocities agree well to approximately the 30 percent radius position. 
Beyond the 30 percent radius position, the computer program overpredicts 
the axial velocity. The computer program modeled the wing vortex system, 
but the actual wing surface geometry was not paneled. Thus, wing thick- 
ness effects were not introduced into the calculations. Any retardation 
of the axial velocity along the horizontal centerline which may exist 
due to wing thickness blockage effects would not be reflected in the 
computer predictions. This may account for the overpredictions of axial 
velocity in Figure 21. 

It is apparent from Figures 18 through 21 that the potential flow 
computer program satisfactorily calculated the flow field at the propeller 
plane of the twin-engine aircraft despite the necessary deletion of a 
large portion of the airframe geometry from the computer model. 

4.3 Nacelle Inlet Inflow Velocity Effects on the Flow Field at the Left 

Propeller Plane 

Additional computer predictions at the twin-engine aircraft propeller 
plane were made in which the nacelle inlet velocity ratio, F/V, was 
changed from a value of 0.29, used in the wind tunnel tests, to a value 
of 0.10. Figure 22 demonstrates the effects of varying the inflow velo- 
city ratio. Plotted are computed radial distributions of angle of rota- 
tional flow along the inboard horizontal centerline (90 degree azimuth) 
of the left propeller plane. As the figure shows, changing F/V produces 
large changes in the flow field only at the hub region directly upstream 
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of the nacelle Inlet, However, the flow field is nearly unaltered at 
radii greater than the inlet radius. These results match the experimental 
findings of Roberts and Yaggy (15). 

That changing of the value of F/V produces large flow field changes 
only in the hub region, where propeller blade loadings are minimal, is 
fortunate. This result indicates that the actual choice of F/V value 
used in the computer simulation of inlet inflow should not be critical 
in regard to its ultimate effect on calculated blade loadings, because 
changing F/V will provide major flow prediction changes only near the 
hub such that changes in loading and vibration of a blade operating in 
the flow field will be slight. 

4.4 Geometry- Induced Contributions Versus Propeller Plane Inclination- 

Induced Contributions to the Flew Field at the Left Propeller Plane 

The flow field at the propeller plane is composed of two contribu- 
tions. These are the contribution due to pure propeller plane inclination 
and that, due to the influence of the airframe. It is of interest to 
examine these two contributions to gain Insights into the relative sig- 
nifigance of each. 

At the plane of an isolated propeller inclined to an angle of attack, 
cig, from the free stream, the free stream velocity contributes azimuthally 
varying radial and tangential velocity components as shown in Figures 23a 
and 23b. An expression for the azimuthal variation of flow angularity, 
in terms of angle of rotational flow 6, at the isolated inclined propeller 
plane may be obtained by inserting the velocity components of Figures 23a 
and 23b. into Equation (48)'. The resulting expression applicable to an 
isolated inclined propeller plane is 

0 - tan (-tan a G sin ij>) 


(54) 
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Flow field predictions produced by the computer program contain both 

| 

the air frame- induced and propeller plane inclination-induced contributions. ] 

i 

A 

By comparing the flow angularity distribution at an isolated inclined j! 

jj 

propeller plane, Equation (54), with the corresponding distribution ob- ij 

l\ 

a 

tained from the computer program predictions incorporating both the in- 
clination and airframe effects, the airframe-induced flow contribution 

\i 

u 

contained in the computer results is highlighted. ji 

Figure 24 presents such a comparison applied to the left propeller I! 

SI 

plane of the twin-engine airplane. In this figure the azimuthal varia- 

il 

tions in angle of rotational flow at the 75 percent radial station are 

| 

presented for both high and low airct'xft angles of attack. As the figure 

If 

indicates, for an aircraft angle of attack of two degrees, the propeller 

plane inclination, a G , is zero. Thus, the isolated propeller plane, in 

the absence of airframe effects, experiences no flow angularity. The 

corresponding curve produced by the computer program shows that the intro- > 

duction of the wing-fuselage nose-left nacelle geometry at two degrees of 

angle of attack produces a nonsymmetric distribution of flow angularity 

with a peak magnitude of nearly two degrees occurring at the 90 degree 

azimuth position. 

For an aircraft angle of attack of 12 degrees, the isolated propeller 
plane inclination is ten degrees. With a ten degree inclination, flow 
angularity for the isolated propeller plane, as shown in Figure 24, is 
symmetrically distributed in a periodic fashion about the azimuth, with a 
maximum flow angularity magnitude of ten degrees at the horizontal center- 
line azimuths of 90 and 2/0 degrees. The corresponding curve calculated 
by the computer program indicates that the addition of the airframe geo- 
metry at an angle of attack of 12 degrees produces greater and more dis- 
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torted flow angularity. Geometry influences have shifted the flow dis- 
tribution curve to the left. Increased upwash introduced by the airframe 
geometry is reflected by the greater amplitudes of 9 in the computer pro- 
gram prediction. Also, the airframe induces more upwash over the inboard 
half of the left propeller plane than over the outboard half. This is 
indicated by greater 0 amplitudes in the vicinity of the 90 degree azimuth 
than in the vicinity of the 270 degree azimuth. 

Figure 25, like Figure 24, also compares the effects of geometry 
influences and propeller plane Inclination Influences. In this figure, 
radial distributions of angle of rotational flow along the Inboard hori- 
zontal centerline, 90 degree azimuth, on the left propeller plane are 
plotted. Because the isolated inclined propeller plane experiences no 
radial variation in the flow field, as indicated by Equation (54), the 
flow angularity due to inclination is a constant for each angle of attack 
in Figure 25. For an aircraft angle of attack of two degrees, the pro- 
peller plane inclination with respect to the free stream is zero, and the 
isolated propeller plane experiences no flow angularity, as indicated by 
the uppermost plot in Figure 25. As the corresponding computer-predicted 
plot shows, the aircraft geometry at a two degree angle of attack Induces 
upwash, particularly strong near the hub due to the nacelle, which pro- 
duces peak flow angularity magnitudes of nearly 5,5 degrees. 

As shown in Figure 25, for an aircraft angle of attack of 12 degrees, 
the isolated propeller plane inclination is ten degrees. Thus, a blade 
in this isolated propeller plane would experience a ten degree flow angu- 
larity at the horizontal centerline equivalent to a uniform blade pitch- 
increase of ten degrees. The corresponding computer program predictions 
for a 12 degree angle of attack indicate a large increase in upwash due 
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to the presence of the airframe geometry. Increased upwash corresponds 
to negative values of 0 in Figure 25. The nacelle is seen to induce the 
greatest upwash in the hub region, while a less but nearly constant flow 
angularity due mostly to wing-induced upwash is seen at the midrange 
radii. Finally, increased flow angularity is induced by the fuselage 
nose near the tip radius in the lowermost plot in. Figure 25. 

As Figures 24 and 25 demonstrate, the influence of the airframe 
geometry on the flow at the propeller plane is very evident. The inclu- 
sion of geometry effects yields more flow angularity and flow field dis- 
tortion than exists when only the propeller plane inclination-induced 
flow is present. 


r 
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Chapter 5 

COMPUTER-PREDICTED FLOW FIELD RESULTS AT THE PROPELLER PLANE 
OF THE PIPER CHEROKEE PA-28-180 AIRPLANE 

5.1 Cherokee 180 Airplane Geometry, Body Paneling, and Run Conditions 

A Piper Cherokee PA-28-180 is owned by the Aerospace Engineering 
Department of The Pennsylvania State University. Because this airplane 
was readily available for research use, computer predictions of the flow 
at the propeller plane (propeller removed) were made for it. 

In this chapter, results of parametric studies are presented which 
show the effects of cowl inlet inflow, wing lift, wing dihedral, and aft 
fuselage geometry deletion on the computed potential flow field at the 
propeller plane. Additionally, a series of figures is included which 
provides a mapping of the computed flow field over the entire propeller 
plane of the baseline Cherokee 180 configuration for a range of fuselage 
angles of attack. 

Figure 26 presents the geometric characteristics of the single- 
engine, fixed-gear Cherokee 180 airplane. Though not indicated in Figure 
26, the propeller is mounted on the airplane such that the propeller 
plane is inclined, with respect to the fuselage, downward and to the 
right by the angular amounts given in Table 2. As noted in Figure 26, 
due to the sideward inclination of the propeller plane, the propeller 
hub position is shifted a distance of 0.0228 propeller radii to the right 
of the fuselage plane of symmetry. 

Computer panel input for all portions of the fuselage aft of the 
firewall was obtained from detailed cross section geometry descriptions 
found on aircraft drawings supplied by the aircraft manufacturer. How- 
ever, all of the cowling surfaces are constructed of fiberglass. As a result, 
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the manufacturer has n6 detailed cross sectional drawings for any of the 
cowling geometry foreward of the firewall. Therefore, in order to gen- 
erate the computer paneling model of the cowling and spinner, it was 
necessary to physically measure the cowling and spinner surface geometry 
of the Univeristy-owned Cherokee airplane. 

Measurement of the cowling and spinner surface coordinates was per- 
formed by triangulation using surveying transits. Three-view drawings 
of the cowling and spinner were made from the measurements. From these 
drawings, the computer paneling model of the cowling and spinner was 
generated. The process of measuring and mapping this geometry is des- 
cribed further in Appendix E. 

None of the tail or wing surfaces were paneled. However, the wing 
was simulated by the horseshoe vortex model used in the computer program. 
Also, none of the landing gear were included in the computer input model. 

Figure 27 presents flight test-measured aircraft lift coefficients 
as a function of fuselage angle of attack. values from this figure 
were assumed to be equivalent to wing lift coefficients and were used as 
program input for wing-on computations. 

Two different aircraft body paneling networks have been created. 

The first network includes the cowl-fuselage combination. The second 
network is the same as the first one, except the fuselage geometry aft 
of the firewall has been replaced by a short, streamlined afterbody. 

This short afterbody is illustrated in Figure 26. The. second panel 
network was used for computer predictions at the propeller plane from 
which the effects of ignoring aft fuselage geometry- induced flow could 


be determined. 
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During preliminary testing of a Cherokee 180 wind tunnel model, a 
cowl inlat inflow velocity ratio, F/V, of 0.2 was measured. This value 
of F/V was accepted as the baseline input value for computer runs. 

Initially, the cowl inlet geometry was precisely paneled producing 
a closed concave panel surface covering the inlet. Unfortunately, the 
iterative solution process of the computer program failed to converge 
because of this concave panel surface. To correct this problem, while 
maintaining a simulation of inlet inflow, the concave inlet paneling was 
replaced by a planar membrane of panels placed over the inlet openinp. 
Each of the panels in this membrane was identified as an inlet panel and 
assigned a value inlet inflow velocity ratio, F/V, to model the inlet 
inflow. 

The spinner paneling also created a problem on early computer runs. 
Initially, the spinner geometry was included and paneled as a discrete 
closed body positioned upstream of the cowl inlet panels. In each 
attempt, the presence of the spinner paneling caused program failure 
manifested by divergence of the iterative solution process. As a result, 
the spinner paneling was deleted from the paneling networks, and all 
computer predictions for the Cherokee 180 were made with the spinner off. 

Table 2 describes the various configurations which were used as 
program input. Additionally, Table 3 lists the values of wing lift 
coefficient, taken from Figure 27, which were used as program input. 
Results from computer runs involving all of these configurations are 
presented in the remainder of this chapter. 

5 . 2 Aft Fuselage- Induced Effects on the Flow Field at the Propeller 

Plane of the Cherokee 180 Airplane 

A study was made to determine how much the accuracy of flow predic- 
tions at the propeller plane of the Cherokee 180 airplane was reduced by 
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deletion of the aft fuselage geometry, that geometry aft of the firewall, 
from the computer panel modeling. Acceptability of neglecting the aft 
fuselage geometry would permit the creation of a computer model containing 
fewer panels, thereby decreasing computing time and cost.. This study 
also was. performed to verify the results of aft fuselage paneling deletion 
previously presented for the arbitrary fuselage in Chapter 3. 

Fuselage-off computer flow predictions were made using the Cherokee 
cowl-short afterbody paneling network and the input parameters of Config- 
uration 6 in Table 2. Effectively, the use of the short afterbody elim- 
inates the windshield geometry and thus reduces the frontal area of the 
paneled configuration. The cowl-short afterbody geometry contains fewer 
panels than the cowl-fuselage geometry and is less expensive to run on the 
computer. Additional fuselage-on flow predictions were made using the 
cowl-fuselage paneling network and the input parameters of Configuration 
1 in. Table 2. 

Comparisons of the fuselage-on predictions with fuselage-off pre- 
dictions were made to determine the significance of the aft fuselage 
contribution to the flow at the propeller plane. Figure 28a compares 
fuselage-on and fuselage-off azimuthal distributions of predicted axial 
velocity. Regardless of the body angle of attack, a, the axial velocity 
magnitude is increased at all azimuths due to the reduced flow blockage 
realized by deletion of the aft fuselage. As the figure shows, the 
fuselage-off overprediction of axial velocity is no more than 1.0 to 2.5 
percent of the free stream velocity. 

Removal of the aft fuselage geometry produces slightly increased 
upwash predictions at the propeller plane. This increase is reflected 
in the fuselage-off tangential velocity distributions in Figure 28b. At 
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the lower body angle of attack, a, deletion of the aft fuselage changes 
the tangential velocity field everywhere except near the vertical center- 
line azimuths of zero and 180 degrees. However, the tangential velocity 
distribution for a body angle of attack of ten degrees is less affected 
by aft fuselage removal, since noticeable changes occur only near the 
horizontal centerline azimuth positions. 

Figures 28c and 28d present the flow angularity distributions cor- 
responding to the previously examined axial and tangential velocity 
predictions for body angles of attack, a, of two and ten degrees, re- 
spectively. At the lower angle of attack. Figure 28c, deletion of the 
aft fuselage geometry results in no more than a 0.3 degree change in 
predicted flow angularity at any azimuth. At the higher angle of attack, 
Figure 28d, fuselage-off computations produce no discernable change in 
the flow angularity distribution. 

The upper vertical centerline of the propeller plane is located 
directly upstream of the windshield surface on the aft fuselage. There- 
fore, presumably, the largest fuselage- induced axial velocity contribu- 
tions at the propeller plane would occur along the upper vertical center- 
line. Figure 29 presents fuselage-on-predicted and fuselage-off-predicted 
radial distributions of axial velocity along the upper vertical centerline 
position (the zero degree azimuth) for high and low fuselage angles of 
attack. Even at this azimuth position, where fuselage influences should 
be strongest, no more than a three percent overprediction of the axial velo- 
city magnitude results from deleting the aft fuselage computer paneling. 

For a typical single-engine aircraft geometry such as the Cherokee 
180, greater accuracy in the flow predictions at the propeller plane is 
achieved by including the aft fuselage geometry in the computer model. 
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However, as the results of this subsection indicate, deletion of the 
aft fuselage geometry from the computer model results in a decrease of 
only a few percent in the accuracy of flow field predictions at the 
propeller plane for a typical single-engine airplane. Specifically, flow 
angularity predictions are negligibly affected, particularly at high 
angles of attack, by aft fuselage geometry deletion. Also, axial velocity 
overpredictions of no more than three percent can be expected should the 
computer paneling exclude the aft fuselage geometry. 

The advantage of computer time and cost savings realized by using a 
fuselage-off panel model may offset the disadvantage of slightly reduced 
flow prediction accuracy. However, if increased accuracy in flow pre- 
dictions is desired, the aft fuselage should be included in the computer 
panel taodel. 

5.3 Wing Lift-Induced and Dihedral-Induced Effects on the Flow Field at 

the Propeller Plane of the Cherokee 180 Airplane 

To determine if it is important to include the wing geometry in the 
computer model of a single-engine airplane for predicting flow at the 
propeller plane, two computer runs were made with the Cherokee 18Q. The 
first computer run was made at both, high and low fuselage angles-of attack 
using the complete cowl-fuselage-wing input geometry (spinner removed) . 

The second run was the. same as the first, except that the wing geometry- was 
deleted from the model. The models for these two runs correspond to 
Configurations .1 and 5, respectively, in Table 2. 

Wing lift-induced upwash in the vicinity of the horizontal center- 
line is the most significant wing contribution to the flow at the. propel- 
ler plane which affects cyclic propeller blade loads and vibration. Thus, 
comparisons of wing-off-predicted and wing-on-predicted radial dtstribu- 






62 


tions of tangential velocity and angle of rotational flow near the hori- 
zontal centerline would indicate how much the wing-induced upwash contri- 
butes to the flow field. 

Such comparisons are presented in Figures 30a, 30b, and 30c for 
radial distributions of flow quantities along the 90 degree azimuth posi- 
tion (near the horizontal centerline) of the Cherokee propeller plane 
using results from the two computer runs. It should be noted that because 
the Cherokee propeller in oriented with vertical and sideward inclinations, 
oip and 6p, the 90 degree azimuth position is a few degrees below the 
actual horizontal centerline rather than coincident with it. 

At the 90 degree azimuth, negative tangential velocities, as in 
Figure 30a, are indicative of upwash. In Figure 30a, the wing-induced 
upwash. increases the magnitude of the tangential velocity by nearly the 
same amount at all radii for a given fuselage angle of attack, a. This 
magnitude Increase is on the order of two percent of the free stream 
velocity for a fuselage angle of attack of two degrees, corresponding 
to a cruise flight condition. However, for a fuselage angle of attack 
of ten degrees (a high lift, takeoff flight condition), wing-induced 
increase in the magnitude of the tangential velocity is as much, as 6.5 
percent of the free stream velocity. 

Figure 30b presents the wing-on and wing-off predictions of angle of 
rotational flow associated with the velocities in Figure 30a for a fuselage 
angle of attack of two degrees. The wing is seen to increase the flow 
angularity at all radii. A one to two degree increase is observed at 
radii greater than the cowl inlet radius, the important radii where the 
propeller blade is heavily loaded. For a fuselage angle of attack of 
ten degrees. Figure 30c indicates wing-induced flow angularity increases 
of two to four degrees at radii greater than the inlet radius. 
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It is evident from Figures 30a, 30b, and 30c that the wing of a 
typical single-engine airplane such as the Cherokee 180 makes a sizeable 
contribution to the flow field at the propeller plane even for low fuse- 
lage angles of attack. Clearly, to ensure more accuracy in the flow 
predictions at the propeller plane, particularly for high fuselage angles 
of attack, the wing geometry must be included in the computer input model. 

Having determined the necessity of including the wing geometry in 
the computer model, the effect of wing dihedral on the computer-predicted 
flow at the propeller plane is of interest. Concern about dihedral- 
induced effects on the flow field was raised during construction of a 

wind tunnel model of the Cherokee 180 to be used in experimental phases 

of this research, To simplify model construction, it was desirous to 
build the wing without dihedral, but it was first necessary to determine 
if dihedral removal would have a negligible effect upon the flow field at 
the propeller plane, to be experimentally measured. 

In an effort to determine the size of the wing dihedral-induced 
contribution to the flow at the propeller plane, two computer runs were 
made. The first run employed the Cherokee cowl-fuselage-wing model with 
wing dihedral included (Configuration 1 in Table 2) . The second run was 

made using no dihedral (Configuration 4 in Table 2). 

Results of the two computer runs are compared in Figures 31a and 
31b. These figures present radial distributions of flow angularity along 
the 90 degree azimuth position for two different fuselage angles of 
attack. For a fuselage angle of attack of two degrees, Figure 31a indi- 
cates that there is no discernable dihedral-induced contribution to flow 
angularity except at the innermost radii, which are embedded within the 
spinner on the actual aircraft and are not important. However, for a 
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fuselage angle of attack of ten degrees (Figure 31b), flow angularity is 
increased at all radii if wing dihedral is not modeled. Dihedral-induced 
decreases in flow angularity of as much as two degrees are observed near 
the 75 percent radial position. 

Obviously, wing dihedral can be deleted from the computer model for 
low fuselage angles of attack, and little loss in flow prediction accuracy 
will result. However, for high fuselage angles of attack, the high wing 
lift magnifies the dihedral-induced contribution to the flow at the pro- 
peller plane. Hence, the actual wing dihedral should be included in the 
computer model to provide increased accuracy in the flow predictions for 
higher fuselage angles of attack. Therefore, in general, modeling the 
wing dihedral is advisable. 

5 . 4 Cowl Inlet Inflow Velocity Effects on the Flow Field at the Propeller 

Plane of the Cherokee 180 Airplane 

Investigation of inlet inflow effects on the flow at the propeller 
plane first performed in Chapter 4 for the twin-engine airplane was con- 
tinued for the single-engine Cherokee 180 airplane. The study was done 
to determine the sensitivity of flow at the propeller plane to changes in 
cowl inlet inflow velocity. 

Using the Cherokoe cowl-fuselage-wing computer model, three computer 
runs were made, each using a different value of cowl inlet inflow velo- 
city ratio, F/V, f.or input. The first run used the baseline F/V value of 
0.2 (.Gonf iguration 1 in Table 2) which was previously measured in wind 
tunnel tests of the Cherokee 180. In the second run, F/V was decreased 
to 0.1 (Configuration 2 in Table 2). F/V was increased to 0.4 in the 
third run (Conf iguation 3 in Table 2). 

In Figures 32a, 32b, and 32c, results obtained using the various 
inlet inflow velocity ratios are compared for both low and high fuselage 
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angles of attack. Radial distributions of axial velocity obtained using 
the three inlet inflow velocity ratios are presented in Figure 32a. For 
both fuselage angles of attack^ changing F/V is seen to produce signifi- 
cant changes in axial velocity only at radii less than the local cowl 
inlet radius, positions which are directly upstream of the inlet. By 
reducing F/V from the baseline value to a value of 0.1, the cowling be- 
comes more like a bluff body at the inlet face, and axial flow blockage 

A 

at the inner radii on the propeller plane is increased. Conversely, by 
increasing F/V from the baseline value of 0.2 to a value of 0.4, axial 
flow blockage at the propeller plane is greatly reduced. In fact, axial 
flow blockage appears to decrease almost linearly with increasing F/V 
at the inner radii. At radii beyond the cowl inlet radius, locations in 
which the propeller blades are heavily loaded, changing F/V has little 
effect on the axial velocity. In fact, changing F/V produces no discernible 
change in axial velocity at locations beyond the 43 percent radial position. 

Radial distributions of flow angularity obtained using the various 
inlet inflow velocity ratios are presented in Figures 32b and 32c for 
fuselage angles of attack of two and ten degrees, respectively. Reasons 
for the strange oscillatory behavior of the flow angularity distributions 
at the inner radii seen in these two figures are given later in sub- 
section 5.5.3. In both figures, it is seen that changes in F/V greatly 
alter the flow angularity only at the inner radii directly upstream of 
the cowl inlet. Flow angularity increases as F/V decreases. At radii 
beyond the cowl inlet radius, where the propeller blades are most heavily 
loaded, flow angularity is little affected by changes in F/V. Beyond the 
50 percent radial position, changing F/V does not noticeably change the 
flow angularity in either of the Figures 32b or 32c. 


66 

The major result observed in this study of inflow effects is that 
changing the input value of inflow velocity ratio, F./V, greatly affects 
the predicted flow field only at the region of the propeller plane dir- 
ectly upstream of the inlet. This finding agrees with the computer- 
predicted results presented in Chapter 4 for the twin-engine airplane 
and agrees with, the experimental results of Roberts and Yaggy (15) . 

The flow field at the hub region directly upstream of the cowl inlet 
has much less affect on the overall propeller vibration and performance 
than does the flow field outside of the hub region. Therefore, any 
inaccuracies in the flow field predictions at the hub region are not 
extremely distressing. Thus, insofar as' overall propeller performance 
and vibration are affected by the flow at the propeller plane, the choice 
of inlet inflow velocity ratio, F/V, used to make the flow predictions 
is not extremely critical. 

To maximize flow prediction accuracy, the actual value of inlet 
Inflow velocity should be specified in the computer program input. But, 
if the actual value of F/V is not known, an approximation of the value 
should suffice. 

5.5 Baseline Flow Predictions Over the Entire Propeller Plane of the 

Cherokee 180 Airplane 

5.5.1 Model used, purpose of presenting the predictions, and organization 
of the data 

In this final section, this chapter presents a series of figures 
providing a complete description of the flow field over the entire pro- 
peller plane of the Piper Cherokee 180 based on computer predictions. 

The predictions were made using the baseline cowl -wing-fuselage (spinner- 
off) computer model (Configuration 1 in Table 2) . 
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Effects of aft fuselage modeling, wing modeling, and cowl inlet 
inflow velocity on the flow predictions at the propeller plane of the 
Cherokee 180, a typical single-engine airplane, were examined in pre- 
vious sections of this chapter. To obtain the most accurate flow pre- 
dictions, it was shown that the computer model must; include the 
aft fuselage paneling; include the wing (and dihedral); and employ 
a value of inlet inflow velocity closely matching the expected in-flight 
value. 

Because the baseline Cherokee computer model incorporates the afore- 
mentioned characteristics, it is the most realistic of the six model 
configurations listed in Table 2. Thus, of all flow predictions from 
the six models, those made using the baseline model are the most accurate. 
Hence, the flow field predictions obtained using the baseline model were 
selected for presentation in this section. 

There is a two-fold purpose in presenting these computer predictions. 
First, the information is presented as a general example, of the flow field 
at the propeller plane of an actual single-engine light aircraft. Second, 
although there are no experimental data currently available for comparison 
with, these predictions to guage the computer prediction accuracy, these 
predictions are also presented to serve as a data base for use during 
subsequent experimental phases of the research project of which the work 
described in this thesis is but a part. 

The series of flow predictions has been divided into two groups 
presented in the following two subsections, the first including only the 
azimuthal distributions of flow parameters, and the second including only 
the. radial distributions. 
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5.5.2 Computed azimuthal distributions 

Figures 33, 34, 35, and 36 present azimuthal distributions of flow 
field parameters at respective radial positions, r/R, of 0.35, 0.5, 0.75, 
and 1,0 on the propeller plane of the Cherokee 180. Parts a, b, and c 
of each of the four figures consist of distributions of axial velocity, 
tangential velocity, and angle of rotational flow, respectively. 

It is instructive to examine these predicted distributions by dis- 
cussing a single flow parameter at a time, plotted on the same part of 
each of the four figures. In this manner, all axial velocity distribu- 
tions will be examined, then all the tangential velocity distributions, 
and last, all the distributions of angle of rotational flow. 

First, consider the axial velocity distributions. Figures 33a, 34a, 
35a, and 36a. At any fixed radius, the azimuthal variation of axial 
velocity magnitude is seen to increase with increasing fuselage angle 
of attack. However, the shapes of the velocity distributions at any par- 
ticular radius are similar for all fuselage angles of attack. With in- 
creasing radius, the axial velocity distributions flatten and deviate 
less from the freestream velocity; thus, a reduction in body-induced flow 
blockage occurs as radius increases. 

For any fixed radius, Figures 33a, 34a, 35a, or 36a, minimum axial 
flow blockage occurs in the region along the upper vertical centerline 

rd? 

(zero degree azimuth position) of the propeller plane for any fuselage 
angle of attack. This region extends up and away from the cowling. In 
fact, at higher fuselage angles of attack, some axial flow acceleration 
occurs at the upper vertical centerline region. This accelerated flow 
region encompasses a larger azimuth range as radius increases. 
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Regardless of radial position, axial flow blockage increases at 

azimuths away from the upper vertical centerline region. In Figures 34a, 

35a, and 36a, greatest flow blockage occurs in a region centered about 

the lower vertical centerline (180 degree azimuth position). At the 0.35 

radial station, Figure 33a, two flow blockage peaks occur, one near the 

120 degree azimuth position and the other near the 240 degree azimuth 

position. At this inner radius, these two azimuth positions are directly 

upstream of and nearest to the cowl surface, hence the peak flow blockage 

at these locations. Note, however, that due to the installed angle of 

sidetilt, 8 , of the propeller plane to the right, the left half-plane 
P 

is positioned slightly farther from the cowl than is the right half-plane. 
Thus, in Figure 33a, the flow blockage near the 120 degree azimuth (on 
the right half-plane) is greater than that near the 240 degree azimuth. 
Similarly, for all the distributions in Figures 34a, 35a, and 36a, maxi- 
mum axial flow blockage occurs at azimuths on the right half-plane due 
to the installed angle of sidetilt, 8 p . 

Second, consider the tangential velocity distributions. Figures 33b, 
34b, 35b, and 36b. Ideally, as previously seen in Figure 23b, for an 
isolated propeller plane with an inclination, a^, from *:he free stream, 
the tangential velocity varies sinusoidally with azimuth and does not vary 
with radius. However, the tangential velocity distributions for the 
Cherokee propeller plane are distorted sinusoids which also vary with 
radius . 

At the 35 percent radial position, Figure 33b, small oscillatory 
irregularities exist on each of the tangential velocity plots. These 
irregularities are primarily a result of cowl-induced effects, which 
quickly diminish with increasing distance from the cowl. That these 

c - V. 
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small Irregularities are cowl-induced is evident, since these irregu- 
larities are not observed in the tangential velocity plots for the outer 
radii, Figures 34b, 35b, and 36b, 

At the upper and lower vertical centerline positions (zero and 180 

degree azimuths, respectively), the tangential velocities would be zero 

if the hub of the propeller plane was located on the fuselage plane of 

symmetry and if there was no installed sidetilt, gp, of the propeller 

plane. However, the hub of the Cherokee propeller plane is offset to the 

right of the plane of symmetry due to the installed sidetilt, g p , of the 

propeller plane. As a result of the offset, the upper and lower vertical 

centerline positions are immersed in a cowl-induced right-directed side- 

wash. Additionally, a right-directed sidewash only due to the right 

sideward inclination, g , of the propeller plane exists at the vertical 

P 

centerline positions. The net sidewash. corresponds to positive and nega- 
tive tangential velocity components at the zero degree and 180 degree 
azimuths, respectively, noted in Figures 33b, 34b, 35b, and 36b, 

Peaks in the tangential velocity profiles at the 35 percent radial 
position, Figure 33b, are much steeper and sharper than the peaks in the 
profiles at the other radii. At the 35 percent radial station, narrow 
regions centered about the 75 and 285 degree azimuths, where the velocity 
peaks occur, lie directly upstream of and are very close to the inlet 
face of the cowl. Thus, cowl-induced effects are intense in the two 
narrow regions. Hence, the peaks in the velocity profiles in Figure 33b 
tend to be very sharp. However, at the 50, 75, and 100 percent radial 
positions, there are no regions on the propeller plane which lie. directly 
upstream of or as close to the inlet face of the cowl. Cowl-induced ef- 
fects are small or negligible at these radii. Hence, the peaks in the 
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velocity profiles in Figures 34b., 35b, and 36b assume rounder, more 
sinusoidal shapes. 

For a fixed fuselage angle of attack, small shifts in the azimuthal 
location of the tangential velocity peaks in Figures 33b, 34b, 35b, and 
36b occur as radial position changes. Partly because of wing dihedral, 
the azimuthal position of the maximum wing- induced flow changes as radius 
changes. This dihedral effect, combined with the effect of decreasing 
cowl-induced flow at larger radii, produces the azimuthal shifts in the 
locations of the tangential velocity peaks as radius changes. 

Also, for any fixed radius, shifts in the azimuthal location of the 
tangential velocity peaks in Figures 33b, 34b, 35b, and 36b occur as 
fuselage angle of attack varies. Apparently, changes in wing lift as 
well as changes in cowl-induced flow with changing fuselage angle of 
attack also contribute to azimuthal shifts in the locations of tangential 
velocity peaks. 

Upwash at an isolated inclined propeller plane is greatest along 
the right and left horizontal centerlines (9.0 and 270 degree azimuths, 
respectively). Thus, as Figure 23b indicates, tangential velocity is 
minimum at the 90 degree azimuth and maximum at the 270 degree azimuth, 
for an isolated propeller plane. However, at the propeller plane of the 
Cherokee 180, Figures 33b, 34b, 35b, and 36b, the peak minimum and peak 
maximum tangential velocities are not azimuthally positioned 18Q degrees 
apart and occur near but not on the horizontal centerline. This is pri- 
marily due to the position of the hub of the propeller plane relative to 
the cowling. 

Additionally, the peak magnitudes at the left half-plane and right 
half-plane azimuth positions in each of the plots in Figures 33b, 34b, 
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35b, and 36b are not equal. This inequality is due to the installed 
propeller plane sidetilt angle, B , Due to B p . the region near the right 
horizontal centerline is slightly closer to the airframe than is the region 
near the left horizontal centerline. Upwash tends to be greater on the 
right half-plane than on the left. Consequently, in each of the plotted 
distributions, the peak positive tangential velocity is slightly larger 
in magnitude than is the peak negative velocity. 

Comparing Figures 33b, 34b, 35b, and 36b, it is seen that the cowl- 
induced and wing-induced effects on the tangential velocity distributions 
for a fiaed fuselage angle of attack diminish with increasing radius. At 
the 35 percent radial station, Figure 33b, wing-induced and cowl-induced 
upwash are greatest and create the largest peak magnitudes of tangential 
velocity, As radius increases (Figures 34b, 35b, and 36b, respectively), 
peak tangential velocity magnitudes decrease. The decrease is mainly due 
to decreased cowl-induced upwash, though it also is due to the slight 
decrease in wing-induced upwash at the outer radii. However, wing-induced 
upwash remains a significant contributor to the tangential velocity varia- 
tion even at the tip radius (Figure 36b) . In Figure 36b, the wing-induced 
effects are clearly evident in the tangential velocity distribution for 
the fuselage angle of attack of four degrees. At this fuselage angle of 
attack, the propeller plane itself has no vertical inclination from the 
free stream, so the tangential velocity profile in Figure 36b for this 
fuselage angle of attack contains no vertical inclination-induced contri- 
butions. Hence, the sinusoidal variation of this tangential velocity 
profile is due almost entirely to wing-induced upwash at the 100 percent 
radial position. 
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With the presence of the airframe-induced upwash, the peak magni- 
tudes of tangential velocity at the propeller plane of the Cherokee 180 
are significantly greater than the peak magnitudes which would exist at 
the same propeller plane operating in isolation, out of the influences 
of the airframe. The air frame- induced increases in tangential velocity 
can be demonstrated by examples in which a comparison is made between the 
tangential velocities at the propeller plane of the Cherokee 180, Figures 
33b, 34b, 35b, and 36b and the tangential velocities at an isolated 
propeller plane having the same inclination from the free stream. 

Take for one example an isolated propeller plane having a geometric 
angle of attack, ctg, of six degrees from the free stream. Because ct^ is 
-4 degrees and 0^ is three degrees for the Cherokee propeller plane, Og 
for the Cherokee propeller plane is approximately six degrees when the 
fuselage angle of attack, a, is ten degrees. Using the expression for 
v t p in Figure 23b, the minimum and maximum tangential velocity ratios for 
the isolated propeller plane with oig equal to six degrees are -0.1455 and 
+0.1455, respectively. However, as seen in Figure 33b, 34b, 35b, and 36b 
for a fuselage angle of attack of ten degrees, the corresponding Cherokee 
propeller plane experiences peak magnitudes of tangential velocity ratio 
far in excess of 0.1455 at all radii. The excesses are due entirely to 
airframe- induced effects. 

For a second example, consider an isolated inclined propeller plane 
which is perpendicular to the free stream. That is, ctg is zero degrees. 
For this isolated propeller plane, the .tangential velocity is everywhere 
equal to zero. The corresponding Cherokee propeller plane is the one for 
which the fuselage angle of attack, a, is four degrees, In contrast to 
the nonexistent tangential velocities on the isolated propeller plane. 
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the plots in Figures 33b, 34b, 35b, and 36b for a equal to four degrees 
indicate nonzero tangential velocities on the propeller plane of the 
Cherokee. Though sidetilt, 0^, contributes- somewhat to these nonzero 
tangential velocities, most of the nonzero contributions are due to the 
airframe-induced flow. 

Third and finally in this subsection, consider the azimuthal distri- 
butions of angle of rotational flow, 6, in Figures 33c, 34c, 35c, and 
36c. These figures present azimuthal variations of flow angularity cor- 
responding to the previously examined axial and tangential velocity 
variations. 

As a reminder, 0 is a function of the ratio of local tangential to 
local axial velocity as given in Equation (48) in Chapter 2 and is an 
angle lying in the plane of the local section of a propeller blade rota- 
ting in the propeller plane. For a given propeller rotational speed, 
the local blade section angle of attack decreases as 0 increases as 
Figure 8 indicates. Thus, Figures 33c, 34c, 35c, and 36c each qualita- 
tively indicate the azimuthal variation in angle of attack of blade sec- 
tions at the given radius on the Cherokee propeller operating in the 
flowfield for a given fuselage angle of attack. 

As indicated in Figures 33c, 34c, 35c, and 36c for any fixed radius, 
azimuthal variation of flow angularity increases with, increasing fuselage 
angle of attack. These increases in azimuthal variation of flow angu- 
larity are due to the combined effects of increased airframe-induced 
upwash and increased propeller plane inclination resulting from increased 
fuselage angle of attack. 

Comparing Figures 33c, 34c, 35c, and 36c, the azimuthal variation 
of flow angularity increases as radial position decreases. Increasing 
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cowl-induced upwash with decreasing radius is the primary cause of this 
trend in the flow angularity distributions. Thus, in the plotted distri- 
butions, a peak 6 magnitude as large as 35 degrees is observed in the 
distribution at the 35 percent radial position for a fuselage angle of 
attack of ten degrees. Conversely, a peak 0 magnitude as small ns three 
degrees is observed at the 100 percent radial position for a fuselage 
angle of attack of two degrees. 

As seen in Figures 33c, 34c, 35c, and 36c, the combined effects of 
cowl-induced flow, wing-induced upwash and propeller plane inclination 
shift the azimuthal positions of the maximum and minimum 0 values on the 
plotted distributions as changes in radial position and fuselage angle 
of attack occur. For an isolated inclined propeller plane, experiencing 
no airframe- induced effects, maximum and minimum 0 values would occur 
exactly at the left and right horizontal centerlines, respectively (the 
270 and 90 degree azimuths, respectively). On all the plotted 0 distri- 
butions at the Cherokee propeller plane with the exception of the plots 
in Figures 35c and 36c for angles of attack of two and four degrees, the 
maximum and minimum 0 values , though shifting slightly in azimuthal posi- 
tion, remain in the vicinity of the left and right horizontal centerline, 
respectively. 

However, in the distributions of 0 for fuselage angles of attack of 
two and four degrees at the two outer radii, Figures 35c and 36c, the 
maximum and minimum values of 0 do not occur near the horizontal center- 
line. Instead, they occur near the upper and lower vertical centerlines, 
respectively (near the zero and 180 degree azimuths, respectively). These 
shifts of maximum and minimum 0 values to the vertical centerline posi- 
tions are due to the effect of the installed sideward inclination, 8 , 
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of the propeller plane. At the outer radii at fuselage angles of attack 
of two and four degrees, the effect of sideward inclination, 8^, of the 
propeller plane outweighs the effects of wing-induced upwash., cowl-induced 
flow, and vertical inclination, a^, of the propeller plane. 

The following example, involving an isolated inclined propeller plane, 
shows the significance of the airframe-induced contribution to the flow 
angularity at the propeller plane of the Cherokee. On an isolated propel- 
ler plane having an inclination, a G , of six degrees from the free stream, 
a peak magnitude of 9 equal to six degrees occurs at the 9,0 and 270 degree 
azimuths for any radius, as computed using Equation ($41 in Chapter 4, 

For a fuselage angle of attack, a, of ten degrees, the propeller plane 
of the Cherokee has an inclination, Oq, from the free stream of approxi- 
mately six degrees. Therefore, 6 values on the Isolated propeller plane 
having an inclination, Oq, of six degrees when compared with the 9 values 
observed on the propeller plane of the Cherokee for a fuselage angle of 
attack, a, of ten degrees will reveal the airframe-induced contribution 
to the observed 9 values on the propeller plane of the Cherokee, For a 
fuselage angle of attack, a, of ten degrees, distributions in Figures 33c, 
34c, 35c, and 36c exhibit peak magnitudes of 9 equal to 34.7, 18.6, 12.7, 
and 10.8 degrees, respectively. Each of these peak magnitudes is greater 
than the six degree value existing on the isolated inclined propeller 
plane solely because the airframe-induced flow at the propeller plane of 
the Cherokee is present. 

A striking illustration of airframe-induced effects on the flow at 
the propeller plane of the Cherokee is obtained if the flow angularity 
increase at the 75 percent radial position, evident in the foregoing exam- 
ple, is reexpressed in terms of an increase in the angle of attack of a 
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propeller blade section operating at that radial position, Operating in 
the propeller plane of the Cherokee, the propeller blade would experience 
a maximum angle of attack when at the 112 degree azimuth, position, since 
the minimum angle of rotational flow, -12,7 degrees-, occurs at that 
azimuth (refer to the distribution in Figure 35c for a fuselage angle 
of attack of ten degrees). However, operating in the isolated inclined 
propeller plane of the foregoing example, the blade section would experi- 
ence a maximum angle of attack when at the 90 degree azimuth, position 
where the minimum angle of rotational flow, —6 degrees, occurs. Assuming 
that in both propeller planes the propeller is operating at an advance 

I 

ratio of 1.2, the decrease in minimum 9 from -6 degrees to -12,7 degrees- 
indicates the following concerning change in blade section angle of at- 
tack: when operating in the flow field at the Cherokee propeller plane, 

Figure 35c, the blade section at the 75 percent radial station attains 
a maximum angle of attack which is 16.96 degrees greater than the maxi- 
mum angle of attack attained when the propeller operates in the isolated 
inclined propeller plane. This increase of 16.96 degrees is entirely 
due to the airframe-induced contribution to the flow angularity at the 
propeller plane of the Cherokee. 

This concludes discussion of the azimuthal distributions of the 
predicted flow parameters at the propeller plane of the Cherokee, In the 
following subsection, the radial distributions of the flow parameters are 
discussed. 

5.5.3 Computed radial distributions 

Figures 37 through 44 present radial distributions of flow- field 
parameters at fixed azimuth positions, t|>, of zero through 315 degrees 
in increments of 45 degrees, respectively, on the propeller plane of the 
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Cherokee 180. Parts a, b, and c of each figure consist of distribution^ 
of axial velocity, tangential velocity, and angle of rotational flow, 
respectively. Predictions for fiv>i values of fuselage angle of attack, 
are included in each figure. 

Indicated in each figure is the local radial position of the peri- 
phery of the cowl inlet. The radius of the inlet varies with azimuth, 
because the inlet opening is elongated, and the hub of the propeller 
plane is positioned to the right of the plane of symmetry. At any azi- 
muth, the local cowl inlet radius is approximately equal to the maximum 
local radial extent of the cowl inlet face, the region on the cowl sur- 
face which is nearest to the propeller plane. Thus, the local radial 
position of the inlet on the figures acts as a reference point, indica- 
ting the radial extent of those cowl surfaces proximal to the propeller 
plane. 

Before discussing in detail each, of the distributions of the various 
flow parameters, the erratic behavior observed at the inner radii on 
almost every plotted distribution in Figure 37a through 44c should be 
discussed. In the hub region at radii less than the local cowl inlet 
radius, the predicted distributions in all figures, with the exception 
of Figure 37a, exhibit unexpected characteristics. The axial yelocity 
profiles have strange peaks or oscillations in this region. The distri- 
butions of tangential velocity and angle of rotational flow are extremely 
oscillatory in the hub region. 

This unexpected erratic nature of the predicted radial distributions 
probably occurs for two reasons. First, for reasons previously given in 
the first section of this chapter, the spinner geometry was omitted from 
the final computer model used to make the flow predictions presented 
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throughout this chapter. Absence of the spinner probably recounts for most 
of the erratic nature of the distributions inside the ten percent radial posi- 
tion. However , on the actual Cherokee 180, the region on the propeller plane 
at positions less than the ten percent radial position is buried within the 
spinner. Thus, the predicted flow inside the ten percent radial position in 
each of Figures 37a through 44c is nonexistant, and in reality should be Ignored. 

Second, possible shortcomings in the paneling mesh at the region of 
the cowl inlet may have caused problems in the predictions at the hub 
region on the propeller plane. The maximum length of some of the panels 
on the inlet face of the cowl in the computer model is greater than or 
nearly equal to the length at (Me spacing between the propeller plane 
and the cowl inlet face. Thus, the paneling mesh at the cowl inlet region 
may be too coarse to produce accurate flow prediction in the hub region 
of the propeller plane. The phenomenon of oscillatory flow distributions 
at the hub region of the propellei plane due to excessive coarseness of 
the paneling mesh was previously observed in the results pertaining to 
flow upstream of a sphere in Chapter 3. The same phenomenon is probably 
occurring in Figure 37a through 44c :> or the Cherokee 18Q, 

Regardless of the causes, any inaccuracies in the flow predictions 
out to a radial position of 20 to 30 percent on the propeller plane of 
the Cherokee are not of extreme concern, because a propeller blade is * 
not heavily loaded at these inner radii. Inaccuracies in the flow pre- 
dictions at these inner radii would have a minimal impact upon loading 
and vibration of propeller blades operating in the flow field presented 
in Figures 37a through 44c. Note, the foregoing discussion on the erratic 
distributions at the inner radii is also applicable to the distributions 
in Figures 30a through. 32c, 
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In the remainder of this subsection, as was done in the previous 
subsection for the azimuthal distributions, the predicted radial distri- 
butions in Figures 32 through 44 will be <;xamined in detail. In turn, 
the group of axial velocity profiles, the group of tangential velocity 
profiles, and the group of profiles of angle of rotational flow will be 
discussed. 

First, consider the radial distributions of axial velocity plotted 
in part a of each of Figures 37 through 44, In each of these figures, 
airframe-induced blockage of axial flow is plainly evident. This block- 
age is greatest at radial positions directly upstream of the cowling, 
radii less than the local cowl ini radius, with axial velocities as low 
as 20 to 35 percent of the free stream velocity existing there. At ail 
azimuths, the airframe-induced blockage of axial flow diminishes rapidly 

as radius increases, since the blockage is nearly entirely due to the 

» 

fuselage (particularly the cowl). Little wing- induced blockage occurs. 

At azimuths of zero, 45, and 315 degrees (Figures 37a, 38a, and 44a, 
respectively) on the upper half-plane, the cowl-induced axial flow block- 
age diminishes more radily than at any other azimuth. Because the hub of 
the Cherokee propeller plane is situated upstream of the upper portion of 
the cowl inlet face, little of the propeller plane at these three azimuths 
lies directly upstream of any cowl surfaces. Hence, cowl-induced blockage 
is only severe at positions less than the 20 to 25 percent radial posi- 
tions in Figures 37a, 38a, and 44a. Beyond the 25 percent radial posi- 
tion, axial velocities in these three figures quickly return, asymptoti- 
cally, to the free stream value. In fact, at higher fuselage angles of 
attack, axial velocity in Figures 37a, 38a, and 44a is accelerated to a 
value greater than the free stream velocity at regions just beyond the 
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cowl inlet radius. However, the axial velocity returns to the free 
stream value as radius continues to increase. 

Because the propeller plane has an installed inclination, 3 , to 

¥ 

the right, the 45 degree azimuth position is rotated toward the cowl and 
wing, while the 315 degree azimuth position is rotated away from the 
cowl and wing. Consequently, the cowl-induced flow blockage at the inner 
radii in Figure 38a is greater than the blockage at the inner radii in 
Figure 44a. Conversely, at the outer radii, flow blockage is less in 
Figure 38a than in Figure 44a. This is due to the existence of slightly 
higher wing-induced axial flow at the 45 degree azimuth than at the 315 
degree azimuth as a consequence of the sideward inclination, 3 p . 

At azimuths of 135, 180, and 225 degrees (Figures 40a, 41a, and 42a, 
respectively) , on the lower half -plane, cowl^induced axial flow blockage 
diminishes much less rapidly with increasing radius than it diminishes 
at azimuths on the upper half-plane, previously examined. Because a 
large area of the lower half-plane lies directly upstream of the cowl, 
significant cowl-induced flow blockage occurs over the entire radius at 
these three azimuths with a three to seven percent reduction of axial 
velocity below the free stream value remaining at the tip radius in each 
of Figures 40a, 4la, and 42a. Though some wing-induced axial flow block- 
age occurs at the outer radii in these three figures, most of the blockage 
at the outer radii is cowl-induced. Also, in contrast to the behavior 
previously noted at azimuths on the upper-half plane, axial flow blockage 
at azimuths of 135, 180, and 225 degrees, on the lower half-plane, in- 
creases with increasing fuselage angle of attack, a. 

Due to the installed propeller plane inclination, 3 p , to the right, 
the 135 degree azimuth position is rotated toward the cowl and wing. 
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while the 225 degree azimuth position is rotated away from the cowl and 
wing. Consequently, axial flow blockage at the 135 degree azimuth posi- 
tion, Figure 40a, is slightly greater than the blockage at the 225 degree 
azimuth position, Figure 42a, at any particular radius for any fuselage 
angle of attack, a. 

At the 90 and 270 degree azimuth positions (near the right and left 
horizontal centerlines, respectively, of the propeller plane) in Figures 
39a and 43a, respectively, cowl-induced axial flow blockage is smaller 
in magnitude and lest extensive (radially) than the blockage at the three 
azimuths on the lower half-plane, previously examined. In Figures 39a 
and 43a, cowl- induced *?,xial flow reduction is significant (more than 
five percent of the free stream velocity) out to approximately the 70 
percent radial position. At the tip radius in Figures 39a and 43a, the 
cowl-induced reduction in axial velocity is not more than two percent of 
the free stream velocity. 

At any particular radius, axial flow blockage at the 90 degree 
azimuth position, Figure 39a, is slightly greater than the blockage at 
the 270 degree azimuth position, Figure 43a. The greater blockage near 
the right horizontal centerline is a consequence of the installed pro- 
peller plane inclination, Bp, which causes the 90 degree azimuth posi- 
tion to lie closer to the cowl than does the 270 degree azimuth position. 

As Figures 39a and 43a also show, the axial velocity profiles near 
the right and left horizontal centerlines of the propeller plane are 
little affected by changes in fuselage angle of attack, a. 

Second, consider the radial distributions of tangential velocity 
plotted in part b of each of Figures 37 through 44. As previously ob- 
served in the radial distributions of axial velocity, airframe-induced 
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tangential velocity contributions also decrease most rapidly with in- 
creasing radius at the zero, 45, and 315 degree azimuths (Figures 37b, 
38b, and 44b, respectively), on the upper half-plane. At azimuths of 
90, 135, 180, 225, and 270 degrees (Figures 39b, 40b, 41b, 42b, and 43b, 
respectively), lying near the horizontal centerline or in the lower half- 
plane, air frame- induced tangential velocity contributions decrease less 
rapidly with increasing radius than at the three azimuths on the upper 
half-plane. At azimuths near the horizontal centerline and on the lower 
half-plane, fairly large airframe-induced tangential velocity contribu- 
tions are observed as far outboard as the 60 percent radial position and 
beyond. Conversely, at azimuths of zero, 45, and 315 degrees, airframe- 
induced contributions (predominantly cowl-induced contributions) are 
fairly large only at positions inboard of the 60 percent radial position. 

Large air frame- induced tangential velocity contributions are more 
extensive on the lower half-plane than on the upper half-plane for two 
reasons. One reason is that most of the frontal area of the cowl lies 
directly downstream of the lower half-plane, so large cowl-induced ef- 
fects persist at outer radii. The other reason is that wing-induced 
upwash, which is strong at all radial positions, heavily contributes to 
the tangential velocties at the 90, 135, 225, and 270 degree azimuth 
locations. 

Ideally, tangential velocity would be nonexistant along the upper 
and lower vertical centerlines (zero and 180 degree azimuths, respective- 
ly) if the hub of the propeller plane was located on the fuselage plane 
of symmetry and if there was no installed sideward inclination, (3^, of 
the propeller plane. However, the hub of the propeller plane of the 
Cherokee is located slightly to the right of the fuselage plane of 
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symmetry, and an inclination, 6 p , to the right does exist. Hence, the 
vertical centerline of the propeller plane lies in a region at which a 
small right-directed sidewash exists. This sldewash field yields small 
positive tangential velocities at the zero degree azimuth position (Fig- 
ure 37b) and yields small negative tangential velocities at the 180 degree 
azimuth position (Figure 41b). In Figure 37b, the tangential velocity 
distribution at the region inboard of the 25 percent radial position is due to 
the cowl-induced sidewash, a consequence of the location of the hub of the 
propeller plane. The nearly constant distribution outboard of the 25 per- 
cent radial position Is largely due to the sidewash component created by 

the installed sidetllt, 8 • In Figure 41b, tangential velocity resulting 

P 

from cowl-induced sidewash persists over nearly the entire radius. The 
nearly constant distribution of velocity noted at the tip radius is a 
result of the installed sidetilt, 8 p » of the propeller plane. 

At azimuth positions of 45, 90, and 135 degrees (Figures 38b, 39b, 
and 40b, respectively), on the right half-plane, negative tangential 
velocities are produced by a positive upwash field. Hence, with increas- 
ing fuselage angles of attack, a, larger negative tangential velocities 
occur at each radius at these three azimuths, because wing-induced and 
propeller plane inclination-induced upwash contributions increase with 
increasing a. Conversely, at azimuth positions of 225, 270, and 315 
degrees (Figures 42b, 43b, and 44b, respectively), on the left half-plane, 
positive tangential velocities are produced by a positive upwash field. 
Hence, tangential velocities increase with increasing fuselage angle of 
attack, a, at each radius at these three azimuths. 

Upwash makes a maximum contribution to tangential velocity along 
the horizontal centerline of the propeller plane, and it makes no con- 
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tribution to tangential velocity along the vertical centerline. Hence, 
as comparison of part b of each of Figures 37 through 44 illustrates, 
the greatest change in tangential velocity per unit change in fuselage 
angle of attack, a, at any particular radius occurs at the 90 and 270 
degree azimuths, which are nearest to the horizontal centerline of the 
propeller plane of the Cherokee. 

Third and finally, consider the radial distributions of angle of 
rotational flow, 0, plotted in part c of each of Figures 37 through 
44. In each of these figures, the flow angularity is greatest in the 
hub region at radii less than the local cowl inlet radius at any azimuth. 
Reasons for the oscillatory, erratic nature of each of the 6 distribu- 
tions in the hub region were previously stated at the beginning of this 
subsection. That negative values of 6 exist on the right half-plane 
(azimuths less than 180 degrees) and positive values of 0 exist on the 
left half-plane indicates the presence of upwash over the entire propel- 
ler plane of the Cherokee at all fuselage angles of attack for which 
flow predictions are presented. 

Comparing all the plots of 0 in part c of each of Figures 37 through 
44, it is seen that flow angularity is generally slightly greater on the 
right half-plane than on the left half-plane. On the right half-plane, 
a magnitude of 0 as 1/Arge as 58.5 degrees is predicted (at the 135 degree 
azimuth position for a fuselage angle of attack, a, of ten degrees). On 
the left half-plane, however, a magnitude of 0 only as large as 53.5 
degrees is predicted (at the 270 degree azimuth position for a fuselage 
angle of attack, a, of ten degrees). That flow angularity is slightly 
greater on the right half-plane than on the left half-plane is attribu- 
table to the existence of the installed right-directed sidetilt, 8 , of 
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the propeller plane of the Cherokee. Due to fT, the right half-plane 
lies slightly closer to the airframe than the left half-plane does. 

Hence, airframe-induced contributions to flow angularity are larger on 
the right half-plane of the propeller than on the left half-plane. 

Figures 37c and Ale clearly show the effect of the installed side- 
tilt angle, 0^, of the propeller plane on the flow angularity at the 
upper and lower vertical centerline azimuths (zero and 180 degrees), 
respectively. At the zero degree azimuth position, cowl-induced flow 
angularity is evident out to the 40 percent radial position and slightly 
varies with changes in fuselage angle of attack, a. This slight varia- 
tion occurs because the hub of the propeller plane is positioned slightly 
to the right of the fuselage plane of symmetry in a region of cowl-induced 
sidewash which slightly changes as fuselage angle of attack, a, changes. 
However, wing-induced flow angularity is negligible at this azimuth. 

Hence, beyond the 40 percent radial position, airframe-induced flow 
angularity is insignificant, and 0 remains fixed at three degrees for 
all fuselage angles of attack due to the installed propeller plane side- 
tilt angle, 0^, of three degrees to the right. Similarly, airframe- 
induced contributions to 0 at the 180 degree azimuth position (Figure 
41c) are negligible beyond the 85 percent radial position. Hence, at 
this radial position, 0 remains fixed at -3 degrees, commensurate with 
the three degree value of propeller plane sidetilt angle, 0^, to the 
right. 

At the six azimuth positions not on the vertical centerline of the 
propeller plane (Figures 38c, 39c, 40c, 42c, 43c, and 44c), flow angu- 
larity increases with increasing fuselage angle of attack, a, at any 
radius. At the inner radii at these azimuths, an increase in fuselage 


87 


angle of attack, ct, of two degrees produces flow angularity increases 
greatly in excess of two degrees, because extra flow contributions in- 
duced by the cowl exist at these radii. At the outer radii, however, a 
two degree increase in fuselage angle of attack, o, produces flow angu- 
larity increases of only 1.5 to 3.5 degrees, because only the wing- 
induced upwash and vertical propeller plane inclination-induced upwash 
are important contributors to flow angularity at these radii. Of course, 
the wing-induced upwash and vertical propeller plane inclination-induced 
upwash make the greatest contribution to angle of rotational flow, 0, at 
the horizontal centerline azimuths. Of all the azimuth positions for 
which results are plotted, consequently, it is at the 90 and 270 degree 
azimuth positions (Figures 39c and 43c) , nearest to the horizontal cen- 
terline of the propeller plane of the Cherokee, that the increase in flow 
angularity per unit increase in fuselage angle of attack, a, at the tip 
radius is greatest. 

The significance of the wing-induced contribution to flow angular- 
ity (0) may be discerned at the tip radius at the 90 and 270 degree 
azimuths (Figures 39c and 43c, respectively) by comparing the flow field 
at the propeller plane of the Cherokee for a given vertical propeller 
plane inclination from the free stream with the flow field at an isolated 
propeller plane having nearly the same vertical inclination from the free 
stream. Due to the combined effects of the installed inclination angles, 
0 ^ and the propeller plane of the Cherokee has a vertical inclination, 
ct„, from the free stream which is approximately four degrees less than 

b 

the corresponding fuselage angle of attack, a. Therefore, the plotted 
flow results at the propeller plane of the Cherokee for a specified value 
of fuselage angle of attack, a, should be compared with flow results at 
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an isolated propeller plane having a vertical inclination, cig, from the 
free stream equal to a — 4 degrees. Then the flow results Jt the outer 
radii of the propeller plane of the Cherokee will differ from the results 
at the isolated propeller plane by an amount only due to the wing-induced 
contribution. At the 90 degree azimuth position on the isolated propel- 
ler plane for vertical inclinations, cu, from the free stream of -2, zero, 
two, four, and six degrees, the angles of rotational flow, 0, are con- 
stant at all radii and are equal to 2.0, 0.0, —2.0, —4.0, and —6.0 degrees, 
respectively (computed using Equation (54) in Chapter 4). However, at 
the tip radius at the 90 degree azimuth position on the propeller plane 
of the Cherokee (Figure 39c) for the corresponding fuselage angles of 
attack, a, of two, four six, eight, and ten degrees, the angles of rota- 
tional flow, 6, are equal to 0.7, -2.0, -4.5, -7.2, and -9.8 degrees, 
respectively. Comparing the results at the two propeller planes, the 
discrepancy between the set of values of 0 for the isolated propeller 
plane and the set of values for the propeller plane of the Cherokee is 
solely due to the wing-induced contribution to the flow angularity at the 
propeller plane of the Cherokee for each fuselage angle of attack, a, A 
similar comparison can be made for the 270 degree azimuth position. At 
the 270 degree azimuth position on the isolated propeller plane for 
vertical inclinations, a G , from the free stream of -2, zero, two, four, 
and six degrees, the angles of rotational flow, 0, are constant at all 
radii and are equal to -2.0, 0.0, 2.0, 4.0 and 6,0 degrees, respectively. 
However, at the tip radius at the 270 degree azimuth on the propeller 
plane of the Cherokee (Figure 43c) for the corresponding fuselage angles 
of attack, a, of two, four, six, eight, and ten degrees, the angles of 
rotational flow, 0, are equal to -0.7, 1.8, 4.6, 7.1, and 9.6 degrees, 
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respectively. Comparing the results at the two propeller planes at the 
270 degree azimuth, the discrepancy between the two sets of values of 0 
is solely due to the wing-induced contribution to the flow angularity at 
the propeller plane of the Cherokee. The foregoing examples for the two 
azimuth positions indicate that the wing-induced contribution to the flow 
angularity is substantial. 

It is apparent from the series of figures presented in this sub- 
section and the preceding subsection that the flow field at the propeller 
plane of the Cherokee 180 airplane is, indeed, very nonuniform. Obviously 
when a propeller is operating in this flow field, the spanwise (radial) 
loading on the blades, is greatly altered, and substantial azimuthal vari- 
ations in blade loading occur, Vibratory stresses in the blades certainly 
are increased because of this nonuniform flow field, .This concludes the 
discussion of the. baseline flow field predictions at the propeller plane 
of the Cherokee 180 airplane. 
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Chapter 6 

CONCLUSIONS AND RECOMMENDATIONS 

A computer program has been developed which computes the three- 
dimensional, steady, incompressible, inviscid, potential flow field at 
a propeller plane (propeller removed) positioned with any installed in- 
clination upstream of an arbitrary airframe geometry. 

Based upon the results of the flow field predictions made for the 
sphere, twin-engine airframe, and single-engine Piper Cherokee airframe 
geometries, four conclusions can be drawn regarding the overall capabili- 
ties of the computer program. Seven more conclusions regarding the 
nature of the flow at the propeller plane and regarding computer program 
input modeling for increased accuracy of the flow predictions at the 
propeller plane can also be made. The latter seven conclusions are based 
upon all the flow prediction results, in general, and upon the results 
of the parametric studies, in particular. 

1. For a simple wingless geometry such as a sphere, the computer 
program yields excellent predictions of the surface flow. Also, accurate 
predictions of the flow field at a propeller plane upstream of such a 
geometry are obtained using the program. 

2. For a typical twin-engine aircraft configuration, the program 
provides very good predictions of velocity profile shapes at the propeller 
plane despite the deletion of aft and remote fuselage and nacelle surfaces 
from the computer model used for the computations. Reasonable predictions 
of velocity magnitude at the propeller plane are obtained using the pro- 
gram for such a configuration based on comparisons with experimental data. 

3. For a typical single-engine aircraft configuration, the program 
provides reasonable predictions of the flow field at the propeller plane. 
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This conclusion, is based upon judgements of the quality of the single- 
engine Cherokee aircraft predictions which are extrapolations from the 
observed quality of the twin-engine aircraft predictions. However, the 
absence of experimental data for comparison precludes making more speci- 
fic conclusions regarding the computer program capabilities for single- 
engine aircraft. 

4. In general, for the typical airframe geometries considered, the 
flow field predictions at any propeller plane obtained from the computer 
program are more accurate at the outer radii than at the inner radii. 
There is room for improvement in the computer predictions at the inner- 
most (hub region) radii. 

5. Flow fields at propeller planes upstream of typical airframe 
geometries do, indeed, exhibit a high degree of flow angularity and 
spatial nonuniformity as the plotted flow predictions for the twin-engine 
and single-engine configurations clearly illustrate. 

6. Though pure inclination of the propeller plane alone is an 
important contributor to the angularity and azimuthal nonuniformity of 
the flow at the propeller plane, the airframe-induced effects are equally 
important or even more important contributors to the flow angularity 

and are much more important contributors to the spatial nonuniformity 
of the flow field. This conclusion is based upon the comparisons between 
the propeller plane flow field results obtained airframe-absent and the 
corresponding results obtained airframe-present for both the single- 
engine and twin-engine aircraft configurations. 

7. The wing must be included in the computer model if greatest 
prediction accuracy from the program is to be achieved. This conclusion 
is based on the computer predictions for the typical aircraft configura- 
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tions which indicate the wing significantly contributes to the flow field 
at the propeller plane, particularly in the region near the horizontal 
centerline. 

8. The actual wing dihedral should always be incorporated in compu- 
ter models to help maximize the accuracy of the flow predictions. This 
conclusion can be Inferred from the observation at the propeller plane 
for higher fuselage angles of attack that flow predictions made using the 
actual wing dihedral appreciably differ from the flow predictions made 
using no dihedral. 

9. For maximum accuracy of the flow field predictions at the pro- 
peller plane, all of the remote component and aft fuselage geometries must 
be included in the paneling input model to the computer program. How- 
ever, fairly good results, which are but a few percent less accurate, can 
be obtained should such remote and aft geometries be excluded from the 
panel model. This is revealed by the studies of aft fuselage effects 
conducted for the simple arbitrary fuselage configuration and for the 
single-engine Cherokee aircraft configuration. The benefits of computing 
cost savings realized by deleting such geometry components from the 
computer model may offset the disadvantages of the slight decrease in 
flow prediction accuracy which results. The user may need to strike a 
compromise between computing costs and prediction accuracy when deciding 
on the size and complexity of paneling models to use. 

10. Given the airframe components to be paneled in the computer 
model, the accuracy of the flow predictions depends on the distribution 
and sizing of surface panels in the model. The paneling mesh must be 
extremely fine and dense over the airframe surfaces immediately aft of 
the propeller plane in order to improve the accuracy of the flow pre- 
dictions at the propeller plane, particularly in the hub region. 
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11. If the true Inlet inflow velocity for program input is not 
known, an approximate value is satisfactory and its use will have only 
a minimal impact on the overall accuracy of the flow predictions at the 
propeller plane. This conclusion is inferred from the observation that 
changes in the specified inlet inflow velocity ratio induce noticeable 
changes in the propeller plane flow field only at the small hub region 
directly upstream of the inlet. Thus, any inaccuracies in the predicted 
flow field due to the input of an approximate value of inlet inflow 
velocity ratio would have only a very small or negligible impact on the 
overall performance of a propeller operating in the flow field. 

Recommendations for future work include the following eight items. 

1. Experimental measurements of the flow field at the propeller 
plane of the Piper Cherokee PA-28-180 aircraft should be obtained and 
compared with baseline flow field predictions from the computer program 
which have been presented in this thesis. Such comparisons would pro- 
vide a definitive check of the program capabilities for typical single- 
engine aircraft configurations. Wind tunnel testing would be the best 
approach, as flow field data at the propeller plane (propeller removed) 
could be gathered and directly compared with the existing baseline pre- 
dictions. Flight testing, however, could not be done propeller-off but 
could be done power-off and would require making measurements at a plane 
just downstream of the actual propeller plane. If flight test data were 
the only kind obtainable, computer predictions for comparison could be 
made at the plane, just downstream of the' actual propeller plane, where 
the flight test data were taken. 

2. The iterative solution method currently employed by the flow 
prediction program. Jumper (33), should be replaced by a direct matrix 
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solution method. Such a program modification would completely eliminate 
solution divergence problems such as those which were encountered during 
this study when attempts were made to model concave inlet surfaces, air-' 
frames including spinners, and remote components of the twin-engine air- 
frame geometry. This program modification could be easily implemented 
by replacing the current subroutines SOLVE and SOLSYM, both of which use 
Gauss-Siedel iteration to solve a system of equations, by new subroutines 
of the same names which perform a direct matrix solution process such as 
matrix inversion or Gaussian elimination. However, as a result of this 
change, a larger computer having more memory would be required to run 
the program. 

3, As the generation of the paneling input and checking for paneling 
errors is the most tedious task faced by a user of the flow prediction 
program, a graphics package capable of generating three-dimensional per- 
spective drawings of the input paneling geometry should be added to the 
current version of the program. Such a graphics package would greatly 
facilitate the process of finding and correcting paneling errors, a pro- 
cess which must be performed before the comparatively costly flow pre- 
diction steps of the computer program can be allowed to proceed. 

4. Flow field predictions for the Piper Cherokee PA-28-180 model 
with the input spinner paneling included should be obtained. Comparisons 
of these predictions with the currently existing baseline predictions 
would indicate if inclusion of the spinner geometry is necessary for 
improved accuracy of flow field predictions in the hub region of the 
propeller plane. With the spinner included as a discrete body in the in- 
put paneling network for the Cherokee 180 airplane, all prediction attempts 
using the current version of the program have, to date, failed due to 
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iterative solution divergence. However, solutions for this airframe- 
spinner combination could be successfully obtained by direct solution 
through the use of the program modified as described in Recommendation 
2, Alternatively, if the airframe-spinngr model were modified by 
fairing the spinner and cowl geometries to form a single body, the cur- 
rent version of the program could be used and Iterative solution diver- 
gence probably would not occur. 

5. Further parametric studies should be done to determine the ef- 
fects on the propeller plane flow field due to changes in the propeller 
plane inclination angles, and 6^; changes in the spacing between the 
cowl and propeller plane or between the nacelle and propeller plane; and 
changes in the lateral and vertical positions of the propeller plane with 
respect to the airframe. Such studies could be done with either the 
twin-engine airplane configuration or the single-engine Cherokee airplane 
configuration. Also, these studies could be done by using either the 
current version of the flow prediction program or a version modified as 
discussed in Recommendation 2. 

6. The current version of the flow prediction program, Jumper (33), 
which solves for the flow field at the propeller plane (propeller re- 
moved) , should be modified to allow the inclusion of propeller interaction 
effects in solving for the airframe-induced flow field at the propeller 
plane. The theory underlying the current flow prediction program permits 
the existence of an onset velocity field, impinging the body, which is 
spatially varying provided the velocity field is steady. By introducing 

a known quasi-steady propeller wake (propeller-induced flow field) as 
part of such a spatially varying onset velocity field impinging the panel 
model of the airframe in the computer program, the resulting flow field 
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predictions at the propeller plane would Incorporate propeller inter- 
action effects. Though the current computer code does not accomodate 
a spatially varying onset velocity field, the program could be changed 
fairly easily to do so by modifying the boundary condition equations 
(generated by subroutines COEFIC and COFSYM) and by modifying the program 
input (subroutine INPUT). This would be the next step toward addressing 
the mutual airframe-propeller interaction problem. 

7. In conjunction with Recommendation 6 for incorporating propeller 
interaction effects in the solution of the airframe-induced flow field 

at the propeller plane, a computer program should be obtained or written 
which will generate a quasi-steady wake downstream of a propeller having 
an arbitrary nonuniform aerodynamic loading. This program would generate 
the spatially nonuniform onset velocity field for use as input to the 
propeller plane flow field prediction program modified as discussed in 
Recommendation 6. Perhaps a vortex lattice method could be employed by 
such a program. 

8. After modifying the flow prediction program to accomodate pro- 
peller interaction effects (Recommendation 6) and after developing a 
propeller wake prediction program (Recommendation 7) , the problem of 
computing the propeller plane flow field, propeller loads, and propeller 
performance for a complete airframe-propeller combination including mut- 
ual propeller-airframe interference effects should be addressed. One 
possibility for addressing this problem is an iterative scheme utilizing 
the propeller performance prediction program, Aljabri (32) ; the flow 
prediction program, Jumper (33), modified to include propeller interac- 
tion effects as described in Recommendation 6: and a propeller wake 
prediction program, Recommendation 7. Such an iterative scheme might 
include the following steps (a through e) . 
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a. For the airframe immersed in a uniform onset velocity 
field (no propeller interaction effects), compute the airframe- 
induced flow field at the propeller plane (propeller removed) 
as has been done in this thesis, but use the modified flow 
prediction program. 

b. Using the resulting nonuniform propeller plane flow 
field from Step a for input, compute the propeller perform- 
ance and loading using the program of Aljabri (32) and compute 
the quasi-stea.ly propeller-induced flow field using a propeller 
wake prediction program. 

c. Create a quasi-steady, spatially nonuniform onset 
velocity field by incorporating the propeller wake results 
from Step b. Then using this nonuniform onset velocity field 
for input to the modified flow prediction program, recompute 
the flow field at the propeller plane. This time, however, 
propeller interaction effects will have been incorporated in 
the flow field solution. 

d. Repeat Step b, but use for input the resulting pro- 
peller plane flow field from Step c. The results would con- 
sist of the computed propeller performance, loading, and wake 
with mutual propeller-airframe interaction effects entering 
into the solutions. Repeat Step c, but use for input the 
results obtained from the repeat of Step b. The result would 
be the computed propeller plane flow field with mutual pro- 
peller-airframe interaction effects entering into the solution. 

e. Repeat, in an iterative fashion, the cycle described 
in Step d until some desired convergence criteria are met. 
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Completion of the iterations would yield the final propeller 
plane flow field; the airframe surface flow field; and the 
propeller performance, loading, and wake results for the air- 
frame-propeller combination with mutual propeller-airframe 
interaction effects accounted for. 

Developing such an iterative procedure would go a long way toward 
achieving the ultimate goal of this ongoing NASA-sponscred project of 
which the work described in this thesis is but a part. That goal is to 
develop the capability of designing aircraft propellers attuned to pos- 
sess optimum performance, vibration, and nclse characteristics when 
operating in combination with a specific airframe configuration. 
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Presented here is the method of partitioning the body surface into 
a network of N panels. Also, formulas for normal unit vector components 
and panel surface areas are presented. These methods and formulas follow 
the techniques used by Woodward, Dvorak, and Geller (27), Additionally, 
formulas for panel control point coordinates are given. All point coor- 
dinates are relative to a single body-fixed Cartesian axis system. 

A. 1 Partitioning of the Body Surface 

A body-fixed Cartesian coordinate system is defined with positive x 
toward the front of the body, positive y to the right, and positive z 
downward. In general, the origin of the coordinate system may be located 
anywhere within or outside the body. However, for a body with left-right 
symmetry, the origin must lie somewhere on the plane of symmetry if the 
advantages of body symmetry are to be obtained. 

The body is divided into a series of cross sections of a constant or 
nearly constant x coordinate. . Cross, sections are specified and numbered 
in sequence from front to aft on the body. 

The Jth cross section is defined by NP discrete periphery points. 

As Figure A. lb shows, for the general nonsymaetric body cross section, 
all NP points are specified in sequence beginning near the top of the 
section and moving clockwise around the section, as viewed toward the 
rear of the body, ending with a repeat of the first point. 

For the special symmetric case. Figure A. la, in which all cross 
sections are symmetric about plane y * 0, only (NP + 1 ) / 2 points need be 
specified in sequence starting at the top centerline and moving around 
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the left (-y) side of the section to the bottom centerline. The computer 
program automatically generates the image points on the other side of the 
plane of symmetry. However, if body sideslip is nonzero, a symmetric 
body must be input in the same manner as a nonsymmetric body. 

Between adjacent cross sections, a ring of panels is generated by 
pairing corresponding points on each section. For example, the first 
and second points on section J are paired with the first and second points 
on section J + 1 to produce the first panel on the ring. Generally, 
panels are four-sided, but triangular panels are created by repeated 
descriptions of a single periphery point. The sequence of specifying 
points is such as to ensure that normal unit vectors, calculated below, 
will be outwardly directed from the body. 

This systematic partitioning produces a network of N panels over 
the entire body (N/2 panels on each side of a symmetric body) . 

A. 2 Panel Normal Unit Vector Components 

Figure A. 2 shows a typical panel created between cross sections J 
and J + 1. Identify the input corner points by 1, 2, 3, and 4 as shown 
on the figure. Define the two diagonal vectors, and B^, given by 

B x - (x 4 - Xj)t + (y A - y x )j + (z^ - z^lc (A.l) 

and 

i 2 “ ( x 3 “ *2^ + ^ y 3 “ y 2^ + ^ z 3 “ z 2^ * ( A - 2 ) 

The vector cross product, (B^ x B 2 ), produces a vector, N, which is 
directed outward from and normal to the panel and is given by 

N-Nt + Nj+h^ , 

x y z 


(A. 3) 
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where the components are 


and 


and where the magnitude Is given by 


Uy 4 - - z 2 ) - 

<y 3 - y 2 >< z 4 " z l )] ’ 

(A. 4a) 

[<x 3 - x 2 )(z a - z x ) “ 

(x^ — x.^) ^ z 3 ” z 2^ * 

(A. 4b) 

[(x 4 - x 1 )(y 3 - y 2 ) - 

(x 3 - x 2 )( y 4 - y x ) ] , 

(A. 4c) 


|n| - (N x 2 + N y 2 + N z V /2 


(A. 5) 


Finally, the outward drawn normal unit vector, n. Is vector N dlvi- 
ded by its own length, j N] . Thus, in terms of the corner points, n is 
given by 

(A. 6) 


n«nt+nj+n^ 
x y z 


with components given by 


and 


N 


N 


N 


N 


N 


m 


(A. 7a) 


(A. 7b) 


(A. 7c) 


where N , N , N , and { N | are given in Equations (A. 4) and (A. 5). 
x y z 


A. 3 Generation of a Flat Quadrilateral Panel 

In general, the four input corner points will not be coplanar. It 
is necessary to have a flat panel for which the surface area is calcula- 
ble, Define the plane of this flat panel as one which is orthogonal to 
unit vector n and contains a point (&, 2) whose coordinates are the 

average of the four input points. Thus, 
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& ■ + x 2 + x 3 + x 4^ * (A. 8a) 

? ■ + y 2 + y 3 + y4> » (a. 8b) 

and 

2 “ -£-(z^ + z 2 z 3 z 4^ * (A. 8c) 

All input corner points 1, 2, 3, and 4 are of equal distance, d, 
from this plane. Distance d is given by 

d » |n x (£ ~ x 2^ + n y (y - y 2 ) + n z ^ z " z 2^l ’ 

where i, and 2 are given in Equation set (A. 8). 

, Input points 1, 2, 3, and 4 are projected distance d onto the new 


plane. Denote the new coplanar point coordinates as xx, yy, and zz. 
Then 


xx k - . 

[ <x lc - \ d) 

9 

k ■ 1 or 4 

(A. 10a) 

, (lt k + n * d > 

9 

k ■ 2 or 3 ; 



f(y k - */) 

t 

k ■ 1 or 4 


yy v ■ • 




(A. 10b) 

K, 

(y k + n y d) 

9 

k ■ 2 or 3 ; 


and 

(z k ' V 0 

9 

k * 1 or 4 


zz, * - 




(A. 10c) 

kC 

(2 k + n z d) 

9 

k ■ 2 or 3 


are the new coplanar points lying on 

the 

corners of the 

flat quadrilateral 


element. Figure A. 3 shows this quadrilateral. 
A. 4 Surface Area of a Quadrilateral P.anel 


Figure A. 4 shows the flat quadrilateral element with the new corner 
points given by Equation set (A. 10). The sides of the quadrilateral have 
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lengths RS, ST, -TU, and UR, and the diagonals have lengths RT and SU. 
Angles at corners 1 and 3 are denoted by ota and S3, respectively. 

In terms of the corner point coordinates, these six lengths are 
given by 

RS - [(xxj^ - xx 3 ) 2 + (yy 1 - yy 3 > 2 + (zz^^ - zz-j) 2 ) 1 ^ 2 , (A.lla) 

ST - [(xxj^ - xx 2 ) 2 + (yy x - yy 2 ) 2 + (zzj^ - zZj,) 2 ]' 1 ^ 2 , (A.llb) 

TU - [(xx^, - xx 2 ) 2 + (yy^ - yy 2 ) 2 + (zz^ - zz 2 ) 2 ] 1/,2 , (A. 11c) 

UR - [(xx 3 - xx 4 ) 2 + (yy 3 - yy^) 2 + (zz 3 - zz^) 2 ] 1 ^ 2 , (A. lid) 

RT - [(xx 2 - xx 3 ) 2 + (yy 2 - yy 3 ) 2 + (zz 2 - zzj) 2 ) 1 ^ 2 , (A. lie) 

and 

SU - [(xx 4 - xx^) 2 + (yy^ - yy x ) 2 + (zz 4 - zz^) 2 ] 1/2 . (A. Ilf) 

Using the cosine law of triangles, the angles are given by 

-1,RS 2 + ST 2 - RT 2 , , k 

““ “ cos 1 2CRSK S T ) 1 (A * 12a) 

and 

-1,DR 2 + RS 2 - su 2 , ,, , 

86 ■ cos [ — HURMRS) — 1 • (A - l2b> 

Finally, the panel area, S, is given by the formula for the area of a 

general plane quadrilateral as follows: 

S - [ (f f - RS) (f f - ST) (f f - TU) (f f - UR) 

- (RS) (ST) (TU) (UR)cos 2 ( aot y - 6 * !] 1 ^ 2 (A.13a) 

where 

ff * y(RS + ST + TU + UR) , (A. 13b) 

and the other quantities are given by Equation sets (A. 11) and (A. 12) . 


A. 5 Control Point on a Quadrilateral Panel 


For the flat quadrilateral panel, the control point (XC, YC, ZC) is 


located at the Intersection of two lines bisecting opposite sides of the 

panel, see Figure A. 5. Define points (G„, G , G ) and (II„, H„, H) as 

x y z x y z 

the midpoints of panel edges RS and TU, respectively. In terms of the 


corner points, 

coordinates of 

points G and H are given by 



G x 


+ XXj) , 

(A. 14a) 


G y 


+ yy 3 ) * 

(A. 14 b) 


G z 

- l (2Z i 

+ zz 3 ) , 

(A. 14c) 


H x 

- 7 <1K 2 

+ ** 4 ) » 

(A.14d) 


H y 

* 7<yy 2 

+ yy 4 ) * 

(A.14e) 

and 






H z 

• J (zz 2 

+ zz^) . 

(A.14f) 

Coordinates of 

control point 

(XC, YC, 

ZC) are then given by 



XC - i(G x 

+ V • 

(A. 15a) 


YC - ±<G y 

+ H y ) , 

(A. 15b) 

and 






ZC - i(G z 

+ v • 

(A. 15c) 


where G , G , G , H , H , and H are given by Equation set (A. 14). 
x y z x y z 

A . 6 Triangular Panels 

As Figures A. 6 a and A. 6 b show, a triangular panel is generated if 
input points 1 and 2 are equal -(this is Case a) or if points 3 and 4' 
are equal (this is Case b) . 
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For triangular panels, the outward normal unit vector, n, is calcu- 
lated as described in Section A. 2. 

The input points themselves are coplanar for triangular panels, so 
the procedure for finding average coplanar corner points, Section A. 3, 
is unnecessary. 

The formulas for calculating surface area and control point coordin- 
ates of triangular panels differ from those used for quadrilateral panels 
and are described in the remaining two subsections. 


A. 6 . 1 Triangular panel surface area 

For each case of a triangular panel, Case a (with the vertex at 
points 1 and 2) and Case b (with the vertex at points 3 and 4), a dif- 
ferent set of expressions for panel area, S, is used. 

First, consider the triangle, Case a, see Figure A. 6a, Define the base 
length as UR and the side lengths as RS and SU. The angle at the vertex 
is aa. In terms of corner point coordinates, the three edge lengths are 
given by 

RS - [(x 3 - x x ) 2 + (y 3 - y x ) 2 + (z 3 - z^ 2 ] 372 , (A.16i) 
SU - [(x 4 - x x ) 2 + (y 4 - y L ) 2 + (z 4 - z^) 2 ] 372 , (A.16b) 

and 

UR - ((x 4 - x 3 ) 2 + (y 4 - y 3 ) 2 + <z 4 - z 3 ) 2 ] 1/2 . (A.16c) 


Using the law of cosines, the angle aa is given by 


aa 


cos 


-1^ (RS) 2 


+ (SU) 2 - (UR) 2 -. 
2(RS) (SU) J 


(A. 17) 


Finally, using the area formula for a general triangle, the surface area, 


S, is 


S - j(RS)(SU)sin(aa) , 


(A. 18) 
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where RS, SU, and aa are given by Equations (A, 16) and (A. 17), respec- 
tively. 

Second, consider the triangle, Case b, see Figure A. 6b, Define the base 
length as ST and the side lengths as RS and TR. The angle at the vertex 

is 08. In terns of corner point coordinates, the three edge lengths 

are given by 

RS - [<x 3 - x x ) 2 + (y 3 - y^ 2 + <z 3 - •, (A. 19a) 

ST - [(x 2 - Xj) 2 + (y 2 - y x ) 2 + (z 2 - z^ 2 ] 1 ^ 2 , (A. 19b) 

and 

TR - [(x 3 - x 2 ) 2 + (y 3 - y 2 ) 2 + (z 3 - z^ 2 ) 1 ^ 2 . (A. 19c) 

Angle 00 is given, using the law of cosines, as 

« ■ 

Then the surface area, S, is given by the area formula for a general 
triangle as follows; 

S - ■|(RS)(TR)sin(00) , (A. 21) 

where RS, TR, and 88 are given, respectively, by Equations (A. 19) and 
(A. 20). 

A. 6. 2 Triangular panel control point 

As with area formulas, the formulas for control point coordinates 
for the triangular panel, Case a, differ slightly from those for the 
other triangular panel, Case b. 

First, consider the triangular panel, Casea, see Figure A. 6a. The 
control point (XC, YC, ZC) is at the area centroid which is located 
midway along line segment DE. Line DE intersects each panel edge one- 
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third the way from the base of the triangle. Thus, points (D^, 

and (E , E , E ) are functions of the corner points as follows: 
x y z 

D x ■ x 3 + 5 <X J - x 3> " 

D y " y 3 + T (y l ‘ y 3 ) • 

D z ' *3 + ' s 3 } - 

E , - x 4 + J (x l ' x 4 ) • 

V ' y 4 + T (y l - y 4> ■ 


D , D ) 
y* i 

(A. 22a) 
(A. 22b) 
(A. 22c) 
(A. 22d) 
(A.22e) 


and 

, E z - z 4 + — z 4 ) » (A. 22f) 

Then in terms of Equation set (A. 22), the coordinates of the control 
point for the triangular panel, Case a, are given by 


XC - y(D x + E x ) , (A. 23a) 

YC - -|(D„ + E ) , (A. 23b) 

+• j y 

and 

ZC “ + E z ) • (A. 23c) 

Second, consider the triangular panel, Caseb, shown in Figure - A. 6b. 

Following the notation of Case a, the control point (XC, YC, ZC) is 

situated at the area centroid, which lies midway along line segment DE. 

Line DE intersects each panel edge one-third the way from the base of 

the triangle. In terms of corner point coordinates, the coordinates of 

points (D , D , D') and (E , E , E ) are given as follows: 
x y 2 * x y z 

D x * x^ + ~-(xj - Xj) , (A. 24a) 

D y ’ y l + 7 (y 3 - y l> • 


(A. 24b) 
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D z " Z 1 + I (z 3 " ®1* * 

E x * x 2 + I (X 3 ‘ x 2> ’ 

E y “ y 2 + I (y 3 " y 2 } ’ 
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(A. 24c) 
(A.24d) 
(A.24e) 


and 

E z “ z 2 + y(z3 - z 2 ) . (A.24f) 

Finally, In terms of Equation set (A. 24), the control point coordinates 
for the triangular panel, Case b, are given by 


and 


* c -i (D * + V • 

YC - |(D y + E y ) , 
ZC-i(D 2 + E 2 ) . 


(A. 25a) 
(A. 25b) 

(A. 25c) 
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Appendix B 

DERIVATION OF FREE STREAM CARTESIAN VELOCITY COMPONENTS 

Denote by V the magnitude of the free stream velocity. Define a 
wind Cartesian axis system, x, y, z, in which V is always directed along 
the negative x axis as in Figure B.l. Next, define by x, y, and z the 
body-fixed Cartesian axes whose origin coincides with the wind axes 
origin. Then for a body having no angle of attack, a, or sideslip, g, f 

the body axes are coincident with the wind axes. This is shown in Figure i 

I 

B.l. | 

For the general case of the body having an angle of attack, a, and 

i! 

sideslip, $, it is necessary to derive expressions for free stream com- 

S 

ponents, u, v, and w, along the x, y, and z axes, respectively. With the 
body and wind axes initially coincident, a pair of rotations is made. 

See Figure B.2. First, a rotation of a about the y axis is made. Denote 
the resulting intermediate set of axes by x*, y*, and z*. Second, a 
rotation of n about the z* axis is made resulting in the body axes,.x, 
y, and z, in their final position with angle of attack, a, and sideslip 
S. Angle n is chosen so as to produce the desired sideslip angle, g, as 

net css 

viewed in the x - y plane. These rotations yield relations between the 
x, y, z system and the x, y, z system. Then derivation with respect to 
time of the axis transformations yields velocities u, v, and w as func- 
tions of V, a, and g. 

Figures B.2a, B.2b, and B.2c show the angular rotations and show the 
resulting projected lengths of all three sets of axes in three different 

I 

views. For scale, the axes have a true arbitrary length, L. | 

It is first necessary to derive the expression for angle n as a j 

i 

U 

ac 9KT 

function of a and g. Refer to Figure B.2a, showing the x - y plane and l 

r; 

\\ 

l \ 

i 
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portraying 3 in its true arc. In this view, axes x*, y*, x, and y pro- 
ject onto an ellipse given by the equation 

[■ ; ■ * - ] 2 + l% 2 - 1 • (B.l) 

L cos ot L 

In Figure B.2a, the projection of axis x ends at point (X,Y) on the 
ellipse. Coordinates of point (X,Y) thus satisfy the following three 


relationships : 


r X ,2 . r Y,2 , 

['■ » 1 " ' J + If] ■ 1 i 

L cos a L 


X ■ L cos a cos ri , 


tan £ ■ ■jj . (B.4) 

Combining Equations (B.2) and (B.3) with Equation (B.4) gives the 
following desired result: 

n ■ tan ^(tan £ cos a) . (B.5) 

As seen in Figure B.2b, the rotation of n about z* yields the fol- 

lowing relations between axis sets x*, y*, z* and x, y, z: 

x* ■ x cos n - y sin q , (B.6a) 

y* ■ x sin n + y cos q , (B.6b) 


Then as shown in Figure B.2c, rotation of a about y yields the 
following relations between axis sets x, y, z and x*, y*, z*: 

x * x* cos a + z* sin a , (1 

y ■ y* t 0 


z * z* cos a - x* sin a . 
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Equation set (B.6) is substituted for x*, y*, and z* in Equation 
set (B.7), and the resulting system is inverted using Cramer’s rule. 
This gives the final relations between the wind axes, x, y, and z, and 
the rotated body axes, x, y, and z, as follows: 



n£ m ■*s 

x * (cos a cos n)x + (sin q)y - (sin a cos n)z » (B.8a) 


MC CCS 

y ■ -(cos a sin n)x + (cos n)y + (sin a 

sin n)z , (B.8b) 

and 

SOS ECS 

z ■ (sin a)x + (cos a)z 

(B.8c) 


The time derivative along each axis gives the velocity component 

parallel to each axis. By definition of the orientation of the wind 

axes 

with the free stream, the following holds: 



« _ v 

dt 

(B.9) 

and 

& . £. . o 
dt dt * 

(B.10) 

Also 

, define the following: 



dx 

u “ dt ’ 

(B.lla) 


v-SX , 
v dt ’ 

(B.llb) 

and 

dz 

w * — . 

dt 

(B. 11c) 


Derivation with respect to time of Equation set 

(B.8) and substitu- 

tion of Equations (B.9), (B.10), and (B.ll) give 



u * -V cos a cos r\ , 

(B. 12a) 


v * V cos a sin n > 

(B.12b) 

and 

w * -V sin a , 

(B. 12c) 
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where n is given by Equation (B.5). 

Equation set (B.12) gives the body Cartesian components of the free 
stream velocity for the body at angle of attack, a, and sideslip, £• 
These components are illustrated in Figure B.3. 
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Appendix C 

COMPUTATION OF VELOCITY COMPONENTS INDUCED BY A STRAIGHT VORTEX 
FILAMENT OF ARBITRARY ORIENTATION 

This appendix derives the formulas used by a computer subroutine to 
calculate the Cartesian velocity components at some specified point, (R, 
2), which are induced by a straight vortex filament of finite length. 
The Biot-Savart law is used to obtain the velocity magnitude, and special 
precautions are used to circumvent numerical difficulties inherent in 
the Biot-Savart law. 

C.l Geometry of the Arbitrary Vortex Filament 

Refer to Figure C.l. The vortex filament of strength, T, starts at 
point 1, having coordinates (k^, z^) , and ends at point 2, having: 

coordinates (x^, y 2 > 2 2 endpoints and T are specified. The vor- 
tex filament coincides with a vector, T, given by 

T - (x 2 - x x )t + (y 2 - y x )j + (z 2 - Zj^k (C.l) 

of length, T, given by 

T - [(x 2 - x x ) 2 + (y 2 - y L ) 2 + - z 1 ) 2 ] 1/2 • (C.2) 

Specify the point, (x, y, z), at which induced velocities are to be 

found. There exists a point C, having coordinates (C , C , C ), coor- 

x y z 

dinates to be determined, positioned on the line containing the vortex 

filament such that a vector, h, from (x, y, z) to (C , C , C ) is per- 

x y z 

pendicular to the filament, h is given by 

h - (C - i)i + (C - y)j + (C - z)ic (C.3) 

x y z 

and has a length, h, given by 

h - [<C x - x) 2 + (C y - h 2 + (C 2 - z) 2 ] 1/2 . 


(C.4) 
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Denote by s a vector between point 1 and point C expressed as 

s - (C x - x^i + (C y - y^l + (C z - z^lc . (C.5) 

Lines D1 and D2 connect point (x, y, z) with the vortex starting point 
and endpoint, respectively. Their lengths are given by the following: 

D1 - [(x - x x ) 2 + (y - y x ) 2 + (z - z 1 ) 2 ] 1/2 (C.6) 


and 

D2 - [ (x - x 2 ) 2 + (y - y 2 ) 2 + (z - z 2 ) 2 ] 1 ^ 2 . (C.7) 

A * 

D1 and D2 intersect the vortex at angles 0^ and 0 2 , respectively. Using 
the law of cosines for triangles, these angles are given in terms of 
previously defined quantities as 


2 2 2 
S T + D1 - D2 
cos 0 i - 2(T)(D1) 


(C.8) 


and 


cos 0 2 *= 


2 2 2 
T + D2 - Dl^ 

2(T) (D2) 


(C. 9) 


It remains to find the coordinates C^, C y , and of point C in 

4 1 

terms of the given geometry. As Figure C.l shows, vectors s and T are 
colinear. Thus, the vector cross product, (s x T), is zero as is each 
component of the vector (s x T). This vector cross product is performed 
using Equations (C.l) and (C.5), and each component of the resulting 
vector is equated with zero. The result is a set of three equations for 
the unknown values of C x » C y , and C^ as follows: 

0 + (z 2 - z x )C y - (y 2 - y x )C z - [y^ - - z^ - y^], (C.10) 

-(z 2 - z x )C x + 0 + (x 2 - x 1 )C z - [z^x, - x 1 ) - (z 2 - Zj)] , (C.l!) 


and 


(y 2 - y^ c x ■ " *l )C y + 0 " ^l ( ^2 ■ ^ ■ ^1 ( *2 ■ *1^' (C - 12 ^ 
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However, the above three equations form an indeterminant system. A 
fourth independent equation is required to solve for C^, C y , and C z » 

^ -fr- 

Because vectors h and T are orthogonal, their scaler product, (h • T) , 
is zero. Using Equations (C.l) and (C.3) in the scaler product and 
equating with zero gives a fourth equation containing unknowns C , C , 

* 7 

and C as follows: 
z 

(x 2 - + (y 2 - y^Cy + (z 2 - z 1 )C z - [x(x 2 - x x ) + y(y 2 - y x ) 

+ z(z 2 - z^)] . (C.13) 


Simultaneous solution of Equations (C.ll), (C.12), and (C.13) for 

unknowns C , C , and C , using Cramer’.s rule, gives the following results: 
x y z 


(x(AA)^ + £.[(BB) Z + (CC)^J * AA[BB(y - y. ) + CC(z « z.))} 

C s r — r , (C.l4a) 

* (AA ) 2 + (BB ) 2 + (CC ) 2 


{yCBBr + y. [ (AA) Z + (CC; ] + BB[AA(x - x.) + CC(z - z.)]> 

C i 5 5 ^-5 » (C.lAb) 

y (AA) + (BB) + (CC) 


and 


(z(CC ) 2 + z. [(AA ) 2 + (BB) 2 ] + CC[AA(x - x. ) + BB(y - y.)]} 

C ± -5 5 , (C.14c) 

Z (AA) Z + (BB) Z + (CC) Z 

where AA - (x 2 - x^) , BB * (y 2 - y.^) , and CC - (z 2 - z^). 


C. 2 Magnitude of the. Total Induced /Velocity 

Denote by W the magnitude of the total velocity induced at point 
(x, y, z) by the vortex filament. When this point lies on the line con- 
taining the filament, length h is zero. In this case, W is fixed at 
zero and the Biot-Savart law is not used. 

When length h is nonzero, the Biot-Savart law is used to calculate 
W. To ensure that W is not unrealistically large, it is checked to see 
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if it is larger than a maximum velocity limit. This velocity limit 
equals 20 percent of the free stream velocity, V. Should W be found 
larger than the limit velocity, W is set 'equal to the limit value, and 
that is the v/’-lue used to obtain the velocity components. 

In summary, velocity W is calculated using one of the following 
equations: 

W «* 0 , if h - 0 ; (C.15a) 

W - ^(cos 0^ + cos 0^) * if h / 0 ; (C.15b) 

or 

W - 0.2V , if W I by Equation (C.15b)] > 0.2V , (C.15b) 
where V is the magnitude of the free stream velocity and all other quan- 
tities have been defined previously. 


C. 3 Induced Velocity Components 

At point (x, y, z) the induced velocity vector, W, is directed 
normal to the plane containing T and h. Define a unit vector, n w , trtiich 
is also normal to the plane containing T and li and has the same direction 

4 " 

as W, see Figure C.l. The vector cross product (h x T) is also a vector 

■4 * 4 “ 4 

normal to the plane containing T and h. Therefore, n can be given by 

w 


± . (h X T) 

n w **• 

|h x T| 


(C.16) 


Substitution of Equations (C.l) and (C.3) into (C.16) yields the follow- 
ing expressions for the components of n^: 

“ Rk [CC(C y ~ ~ BB(C Z ~ 2)] • (C.17a) 

n Wy - r^;[AA(C 2 - z) - CC(C X - £)] , CC .17b) 

°w z * r^[ b B(C x - x) - AACC y - y)] , 


and 


(C.17c) 
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where 

RXR - {[CC(C - y) - BB(C - z)] 2 + [AA(C - z) - CC(C - x)) 2 
y z z X 

+ [BB(C x - £) - AA(C y - y)] 2 } 1/2 , (C.17d) 

and where AA ■ (x 2 - x^), BB ■ (y 2 - y^, CC * (z 2 - z^) , and C x , C y , and 
C are given by Equation set (C.14). Express the induced velocity vector 
in terms of its Cartesian components as follows: 

ft * u^t + VjJ + w . (C.18) 


These components are shown in Figure C.l. Because W coincides with unit 


vector n w , the following holds: 


+ W 

n * — 

w VI 


(C.19) 


Equations (C.15), (C.18), and (C.17) are substituted for W, W, and com- 
ponents of n w , respectively, in Equation (C.19). By comparing like terms 
on each side of the resulting equality, it follows that the desired com- 
ponents of the vortex filament™ induced velocity are given by 


u. 


"Wx w 




(C. 20a) 
(C.20b) 


ana 


w. 


n w z w 


(C.20c) 
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Appendix D 

TRANSFORMATIONS RELATING THE PROPELLER PLANE AXIS SYSTEM 
TO THE BODY-FIXED AXIS SYSTEM 


D.l Coordinate Transformations 

Denote by x, y, and z the Cartesian coordinate axes which remain 
attached to the propeller plane. Positive x is directed along the pro- 
peller axis of rotation in the thrust direction as in Figure D.l. A 
point, P, on the propeller plane is positioned relative to a cylindrical 
coordinate system, r - i^», where azimuth angle, ip, is measured clockwise 
from the negative z axis. Also denote by x 1 , y' , and z' three Cartesian 
axes which remain parallel to the body-fixed axes, x, y, and z, but which 
are displaced to a point (X^ ub> Y^ b> relative to the body-fixed 

axes. Figure D.l shows all the axis systems and shows the propeller 

plane initially without angle of attack, a , or sideslip, B . In this 

P P 

situation the x, y, and z axes are coincident with the x' , y’ , and z' 
axes, respectively. 

For the general case of the propeller plane oriented at some combin- 
ation of angle of attack, a^, and sideslip, Bp, relative to the body, it 
is necessary to express the r - ip coordinate system, attached to the pro- 
peller plane, in terms of the body-fixed axes, x, y, and z. 

With the x, y, and z axes initially coincident with the x' , y\, and 
z’ axes, a pair of rotations is made, see Figure D.2. First, a rotation 

of a about the y* axis is made. Denote the resulting intermediate set 
P 

of axes by x*, y*, and z*. Second, a rotation of n about the z* axis 
is made resulting in the propeller plane-fixed axes, x, y, and z, in 
their final skewed orientation with angle of attack, a^, and sideslip. 

Bp, with respect to the x* , y', and z' axes. Angle Pp is chosen so as 
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to produce the desired sideslip angle, 0^, as viewed in the x' - y‘ plane. 

Once transformations between the skewed x, y, and z axes and the x’, y’, 

and z' axes have been obtained, the final transformations between the 

r - \p and x - y - z axis systems are then obtained. 

Figures D.2a, D.2b, and D.2c show the two angular rotations and show 

the resulting projected lengths of all four sets of Cartesian axes in three 

different views. For scale, the axes have a true arbitrary length, L. 

It is first necessary to derive the expression for angle n as a 

P 

function of and 6^. Refer to Figure D.2a, showing the x' - y* plane 
and portraying 0 p in its true arc. In this view, axes x*, y*, x, and y 
project onto an ellipse given by the equation 

■r^-r-i 2 + - 1 • 

P 

In Figure D.2a, the projection of axis x ends at point (X, Y) on the 
ellipse. Coordinates of point (X, Y) thus satisfy the following three 
relationships : 


[ 


X ,2 + rY.2 

L cos a 1 l L J 


(D.2) 


X ■ L cos a cos ri 
P P 


CD. 3) 


and 

tan 0 * £ . (D.4) 

P x 

Combining Equations (D.2) and (D.3) with Equation (D.4) gives 

n ■ tan” 1- (tan 0 cos a j , (D.5) 

P P P 

which is the desired function relating n to a and 0 . 

P P P 

As seen in Figure D.2b, the rotation of rip about z* yields the fol- 
lowing relationships between axis sets x*, y*, z* and x, y, z: 


x* * x cos n “ y sin n 
p p 


» 


(D.6a) 
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y* ■ x sin n p + Y cos n p , 

(D.6b) 

and 

z* ■ z . 

(D.6c) 


As shown in Figure D.2c, rotation of ct p about y' yields the follow- 

ing relationships between axis sets x' , y* , z 1 and x*, y*, z*: 



x’ ■ x* cos a + z* sin a , 

P P 

(D.7a) 


y* - y* * 

(D.7b) 

and 

z' * z* cos a - x* sin a 
P P 

(D.7c) 


Axes x' , y', and z' are displaced to point Z hub^ 

rela- 

tive 

to the body-fixed axes x, y, and z, so the following applies: 



x ■ x ' + ^Ub • 

(D.8a) 


yy' * \ub • 

(D.8b) 

and 

2 • z ' + z hisb • 

(D,8c) 


Equation set (D.6) is substituted into Equation set (D.7) to elimin- 
ate x*, y*, and z*. Then, substitution into Equation set (D.8) produces 
the following relationships between the body-fixed axes, x, y, and z, 
and the propeller plane-fixed axes, x, y, and z: 

x * (cos a cos n )x - Ceos a sin n )y + (sin a )z + X, , , (P.9a) 

p p p p p nub 

y - (sin n p )x + (cos n p )y + Y hub , (D.9b) 

and 

z * - (sin a cos n )x + (sin a sin n )y + (cos a )z + 7.. . (D.9c) 

p p p p p nub 

As shown in Figure D.3, the cylindrical coordinate system, r - i|i, 

on the propeller plane is related to the propeller plane-fixed x, y, and 



z axes as follows: 
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and 


x » 0 (on the propeller plane) , 
y ■ r sin ^ , 


(D.lOa) 

(D.lOb) 


z * -r cos ^ . GP.lOc) 

Equation set (D.10) is- substituted Into Equation set C D,9. ), Then 
point P, at coordinates (r, ip) on the propeller plane, has corresponding 
Cartesian coordinates (X p , Y p , Z p ) in the body-fixed axis system given 
by 


X 

P 


-r(sin ij/ cos a sin n + cos ^ sin a ) + X^ ub 


Y p - r (sin * cos n p ) + Y hub , 


CD. 11a) 


CD. lib) 


and 


E - r^sin ip f j n o sin n - cos .ip cos o ) + Z. . , CD. 11c) 

p p p p nuD 

where n p is given by Equation (D.5). 

Equation set CD. 11) is the resulting transformation from the propel- 
ler plane-fixed cylindrical coordinates to the body-fixed Cartesian co- 
ordinates for a point on a propeller plane centered at CX^^* ^hub^ 

and tilted with an angle of attack, a p , and sideslip, £ p , relative to the 
body-fixed axes. 


D. 2 Velocity Component Transformations 

It is necessary to derive transformations which, convert Cartesian 
velocity components, u p , v , and w p , at point P on the propeller plane, 
relative to the body-fixed axes, into corresponding cylindrical system 
velocity components, v 3p , v rp , anc * v tp' As shown in Figure D.3, v ap , 
Vj. p , and Vtp are axial, radial, and tangential velocity components-, re- 
spectively, with the axial component directed normal to the propeller 
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plane and positive in the thrust direction (.into the figure) . Tangential 
velocity, v tp> at ^ is in-plane and positive if it follows the right-hand 
rotation of a propeller. v r is in-plane and positive radially outward. 

Equation set (D.9) is inverted using Cramer's rule giving expressions 
relating the propeller plane-fixed axes, x, y, and z to the body-fixed 
axes, x, y, and z, as follows: 

x - (cos a p cos n p ) (x - X hub ) + (sin q p ) (y - * hub ) 

- (sin a cos q ) (z - Z. . ) , (D.12a) 

p p nuD 

7 - -(cos Op sin n p )(x - X^) + (cos n p ) (y - Y hub ) 

+ (sin o sin q ) (z - Z . ) , (D.12b) 

p p hub 

and 

z - (sin o p )(x - X^ uh ) + (cos a p ) (z - Z^) . (D.12c) 

At a point on the propeller nlane, the velocity components, u p , v p , 

and w , relative to the body-fixed Cartesian axis system are given by 
P 

the following time derivatives: 



dx 

u p " dt ’ 

(D.13a) 


dy 

V p “ dt ’ 

(D . 13b) 

and 

dz 

w = -rr . 
p dt 

(D.13c) 

Similarly, 

the velocity components, u p , v p . 

and w , relative to the pro- 
P 


peller plane-fixed Cartesian axis system are given by time derivatives 

(P.14a) 


as follows: 


— dx 
U p " dt 



CD. 14b) 
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and 

w - ~ . (D.l4c) 

p dt 

The time derivative of Equation set (D.12) it taken and Equations (D.13) 
and (D.14) are substituted. Observe that the derivatives of 
and vanish. The result is a set of relations transforming the velo- 
city components, u , v , and w , into components, u , v , and w , as 

P P P P P P 

follows: 

u * (cos o cos n )u + (sin n )v - (sin a cos n )w , (D.15a) 

P P PP PP P PP 

v ■ -(cos a sin n )u + (cos n )v + (sin a sin n )w , (D.15b) 

P PPP PP P PP 

and 


w = (sin a )u + (cos a )w 
P P P P P 


(D.15c) 


As indicated in Figure D.3, propeller plane-fixed Cartesian velocity 
components, u^, v^, and w^, are related to the corresponding axial, radi- 
al, and tangential velocities, v ap* v rp» an< * v tp* at P°^ nt ** as follows: 

v 3 p * ^ , (D.16a) 

v r * (sin \p)v - (cos if>)w , (D.16b) 

PPP 

and 

v t p * (cos 4>) v p + (sin 'l')Wp . (D.16c) 

Finally, Equation set (D.15) is substituted into Equation set (D.16). 

This produces the following relationships: 

Vo * (cos a cos n )u + (sin n )v - (sin a cos n )w , (D.17a) 

P PPP PP PPP 

v r * -(sin a cos \b + cos a sin n. sin t|»)u + (cos n sin ifOv 

r P P P P P P P 

+ (sin a sin n sin ill - cos a cos ip)w , (D.17b) 

P P P P 
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and 

v* m (gin a sin ^ - cos a sin q cos i)/)u + (cos n cos i|<)v 

c p p P P P P P 

+ (cos a p sin ip + sin a p sin cos , (D.17c) 

where ti is given by Equation (D.5). For a point P at position (r, i|>) 
on the propeller plane, Equations (D.17) are the desired relations which 
transform the body-fixed Cartesian velocity components at P, Up, v , and 
Wp, into the corresponding axial, radial, and tangential velocity com- 
ponents, v a , v r , and v t , relative to the propeller plane-fixed cylin- 
P P P 

drical coordinate system. 
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Appendix E 

MEASUREMENT AND MAPPING OF THE PIPER CHEROKEE PA-28-180 COWL 
AND SPINNER SURFACE GEOMETRY 

No manufacturer-supplied detail drawings of the cowl and spinner 
exist for the Piper Cherokee PA-28-180. As these airframe components 
are proximate to the propeller plane, a detailed description of these 
surface shapes was important for creating a good computer paneling model. 
Therefore, to obtain this- description the only recourse was to physically 
measure the cowl and spinner geometry on an actual aircraft. Direct 
measurement of this geometry on The Pennsylvania State University-owned 
Cherokee 180 research airplane, tail number N907PS, was conducted. This 
measurement process is described in this appendix. 

Measurement by triangulation was performed using two precision sur- 
veying transits, both capable of angular measurement to within a toler- 
ance of +20 seconds of arc. Figures' E.la and E.lb schematically illus- 
trate the deployment of the two transits, identified as transit A and 
transit B. 

Parked on a flat surface, the aircraft fuselage was jacked to the 
horizontal position. The horizontal position was determined by placing 
a level on the reference leveling screws located on the exterior of the 
fuselage beneath, the pilot’s window as shown in Figure E.lb, 

Transit A was positioned on the fuselage plane of symmetry at a 
distance of approximately 30 percent wing span forward of the spinner. 

As shown in Figure E.lb, transit A was vertically positioned such that 
the horizontal line of sight intersected the spinner near its tip. This 
intersection point, on the plane of symmetry, was marked and designated 
as the coordinate origin for all measurements. The horizontal distance 
X A between transit A and the origin was precisely measured. 
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Next, using sightings from transit A, transit B was positioned 
directly to the left of the coordinate origin at a precisely measured 
horizontal distance Xg as shown in Figure E.la. Thus, the lines of 
sight of both transits, viewed from above, intersected at a right angle 
at the origin. The height above ground of transit B was slightly less 
than the height of the origin. The vertical distance Ay between the 
origin and transit B was measured. Knowledge of x A » Xg» and Ax completely 
determined the transit positions with respect to the origin, and the 
transits were then ready for use. 

The cowl and spinner surfaces were prepared for measurement by cov- 
ering the left half of both with a mesh of tiny adhesive paper target 
points. These paper targets were visually placed in a series of rows of 
nearly constant x coordinate such that a row of points approximately de- 
fined a body cross section. The target points were closely spaced at 
the forward end of the cowling and at other regions of extreme surface 
curvature. Target points were more thinly spaced at regions on the cowl 
where the surface was less convoluted. 

After placement of the surface targets, actual measurement was begun. 
As shown in Figures E.la and E.lb, a measurement consisted of simultaneous 
sightings on a surface target point by transits A and B. While sighted 
on a target, the horizontal angles and ij/g were recorded. Also, the 
vertical angles and aig were recorded. This process was repeated for 
all target points on the left side of the cowl and spinner. 

Transit B was then moved directly across to the right side of the 
spinner origin point. Additional target points were placed on ducts and 
surface features found only on the right half of the cowl. Further sight- 
ings were made to measure these target points. This concluded the gath~ 
ering of data. 


130 


Using the measured lengths x^> Xg» an d Ax, as well as the four angles 
^A’ ^B' U A* anc * w b’ assoc *- ate( * with a target point, trigonometry was ap- 
plied to calculate the Cartesian coordinates of the target point relative 
to the spinner origin. By using a computer program, this trigonometric 
analysis was rapidly completed for the hundreds of target points. The 
result was a discrete point description of the cowl and spinner geometry 
in Cartesian coordinates. 

Preliminary three-view cowl and spinner drawings were made using the 
measured surface points. As previously stated, the target points had 
been positioned visually in a series of rows, each row lying approximately 
at a constant x coordinate. Of course, analysis of the data revealed 
that points on a given row were not precisely aligned at a single x co- 
ordinate. Therefore, the measured surface points on each row were shifted 
slightly until all rows fell on contours of constant x coordinate and thus 
defined true cross sections. Guided by the preliminary three-view drawings, 
the data point shifting was done carefully using cubic spline interpola- 
tion. 

Final three-view cowl and spinner drawings were made using the shifts 
ed surface points. The front view of these drawings afforded detailed 
descriptions of geometry cross sections. 

Surface coordinates from the final three-view drawings were used to 
generate cowl and spinner input computer paneling. To produce a symmetric 
cowl paneling network, small irregular cowl surface features, such, as 
engine exhaust pipes and small intake ducts, were ignored. 

The manufacturer-defined Cherokee 180 aircraft coordinate origin, 
used in all existing drawings of fuselage geometry, aft of the firewall, 
did not coincide with the arbitrarily chosen cowl-spinner coordinate 
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origin. Therefore, to properly mate the measured cowl-spinner geometry 
with the remaining aircraft geometry, coordinates of all geometry on 
manufacturer-supplied drawings were converted to the measured cowl-spinner 
coordinate system. Thus, all aircraft input geometry, including the cowl, 
fuselage, and spinner paneling, for the flow prediction program was ref- 
erenced to the chosen coordinate origin indicated in Figures E.la and 


E.lb. 
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Table 3 

Piper Cherokee PA-28-180 Aircraft Lift Coefficients Corresponding 
to the Angles of Attack Modeled by the Computer Program 


Fuselage angle of attack, 
a (degrees) 

2.0 

4.0 

6.0 
8.0 

10.0 


C^ (from Figure 27) 

0.385 

0.568 

0.751 

0.934 

1.116 


■l- 
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Figure 3. Front View of a Symmetric Body Panel -Network 
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DISTANCE UPSTREAM OF LEADING EDGE, x/c 


Figure 6. Comparison of "Exact" and Approximate Upwash Distributions 
Upstream of a Flat Plate Airfoil 






poor SSEnv 


MAIN PROGRAM 



CALL COFSYM 


COFSYM 

** 

WINGV 1 


1 



it 

CALL SOLYSM 

1 l-J 

SOLSYM 


VORTEX 

FF 




CALL 

COEFIC 

~ 


COEFIC 

L-^f 


J 

] — 





1 

r 

i 


| CALL 

SOLVE 

— ! 


-U 

SOLVE 

l 

■ 

y — — 

-r 



1 



r 

i 



WINGV 

~W 


VORTEX I 


K - K + lk 1 K - 0 I | ! 

WINGV VORTEX 1 

^ i 

—it . 


CALL VELOCI 





VORTEX 

Z5Z 


tea WINGV 




EULER 


Figure 9. Organization and Structure of the Three-Dimensional 
Potential Flow Computer Program 
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Figure 12. Axial Velocity Distribution Along the Upper Vertical 
Centerline of a Propeller Plane Upstream of a Sphere 
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Figure 14. Axial Velocity Distribution on the Upper Vertical Centerline 
of a Propeller Plane Predicted With a Complete Simple Fuse- 
lage Shape and With the Isolated Cowling 
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Root section NACA 23015 
Tip section NACA 23012 


All dimensions ape in meters 
unless otherwise noted. 


Only the left nacelle and 
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STATION 

NUMBER 

FUSELAGE COORDINATES (METERS) \ 

DISTANCE AFT 

a' 

b' 

c’ 

1 

0.00000 

0.00000 

0.00000 

0.00000 

2 

0.02540 

0.07620 

0.05842 

0.06350 

3 

0.05080 

0.11430 

0.08890 

0.08890 

4 

0.07620 

0.13462 

0.12192 

0.11430 

5 

0.10160 

0.16002 

0.14732 

0.12700 

6 

0.15240 

0.19050 

0.18796 

0.15748 

7 

0.22860 

0. 22352 

0.23876 

0.19304 

8 

0.30480 

0.25400 

0.28702 

0.21590 

9 

0.60960 

0.34290 

0.43180 

0.29210 

10 

0.91440 

0.41656 

0.54610 

0.34544 

11 

1.21920 

0.46228 

0.64770 

0.39370 

12 

1.52400 

0.49530 

0.73152 

0.42164 

13 

1.82880 

0.52832 

0.81280 

0.44704 

14 

2.13360 

0.55372 

0.87376 

0.47244 

15 

2.43840 

0.57404 

0.93472 

0.48768 

16 

2.74320 

0.59182 

0.99060 

0.50292 

17 

2.99720 

0.59944 

1.03378 

0.51054 

18 

3.04800 

0.59944 

1.06680 

0.51054 

19 

3.20040 

0.60706 

1.19380 

0.51308 

20 

3.35280 

0.60960 

1.31064 

0.51562 

21 

3.50520 

0.61214 • 

1.40208 

0.51816 

22 

3.65760 

0.61468 

1.44780 

0.52070 

23 

3.96240 

0.61722 

1.47320 

0.52578 

24 

4.25450 

0.61722 

1.46050 

0.53086 


Figure 17a. Twin-Engine Airplane Fuselage Coordinates for the Portions 
Forward of the Wing Leading Edge (Source: Reference 15) 



152 


ORIGINAL 

0F poor 
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QUALITY 



STATION 

NUMBER 

NACELLE COORDINATES (METERS) ! 

DISTANCE AFT 

r’ 

1 

0.0000 

0.54610 

2 

0.0254 

0.58420 

3 

0.0508 

0.60960 

4 

0.1016 

0.63500 

5 

0.1524 

0.66040 

6 

0.2286 

0.68580 

7 

0.3048 

0.70612 

r 8 

0.4572 

0.72390 

9 

0.6096 

0.73660 

10 

0.7620 

0.74930 

ii 

0.9144 

0.75184 

12 

1.0668 

0.75438 

13 

1.2192 

0.75692 

14 

1.3716 

0,75692 

15 

1,5748 

0.75692 


Figure 17b. Twin-Engine Airplane Nacelle Coordinates for the Portions 
Forward of the Wing Leading Edge (Source: Reference 15) 
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Figure 19. Comparison of 'Computed and Measured Azimuthal Distributions 
of Angle of Rotational Flow at the Left Propeller Plane of 
the Twin-Engine Airplane 
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Figure 20. Comparison of Computed and Measured Radial Distributions of 
Angle of Rotational Flow at the Left Propeller Plane of the 
Twin-Engine Airplane 
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NONDIMENS IONAL RADIAL POSITION, r/R 


Figure 21, Comparison of Computed and Measured Radial Distributions of 
Axial Velocity at the Left Propeller Plane of the Twin- 
Engine Airplane 
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Figure 22. Effect of Change in Nacelle Inlet Inflow Velocity on the 

Radial Distribution of Flow Angularity at the Left Propeller 
Plane of the Twin-Engine Airplane 
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Isolated Inclined Propeller Figure 23b. Velocity Components at an Isolated Inclined 

Plane (Left Side View) Propeller Plane (Forward View). 
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AZIMUTH POSITION, ip (DEGREES) 

Figure 24. Comparison of Azimuthal Distributions of Flow Angularity at 

the Left Propeller Plane of the Twin-Engine Airplane Obtained 
With and Without Airframe Geometry Influences Present 
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Figure 25. Comparison of Radial Distributions of Flow Angularity at the 
Left Propeller Plane of the Twin-Engine Airplane Obtained 
With and Without Airframe Geometry Influences Present 
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Note: Propeller huh position is 0.0228 propeller radii to the right of 

the fuselage plane of symmetry due to propeller plane inclination. 
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Figure 27. Lift Curve for the Piper Cherokee PA-28-180 Airplane 
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Figure 28a. Effect of Neglecting the Piper Cherokee 180 Fuselage Geometry 
(Replaced by a Short Afterbody) on the Computed Azimuthal 
Distribution of Axial Velocity at the Propeller Plane 
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Figure 28b. Effect of Neglecting the Piper Cherokee 180 Fuselage Geometry 
(Replaced by a Short Afterbody), on the Computed Azimuthal 
Distribution of Tangential Velocity at the Propeller Plane 
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Figure 28c, Effect of Neglecting the Piper Cherokee 180 Fuselage Geometry 
(Replaced by a Short Afterbody) on the Computed Azimuthal Distribu- 
tion of Flow Angularity at the Propeller Plane (a * 2° Case) 
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Figure 28d. Effect of Neglecting the Piper Cherokee 180 Fuselage Geometry 
(Replaced by a Short Afterbody) on the Computed Azimuthal Distri- 
bution of Flow Angularity at the Propeller Plane (a « 10® Case) 
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NONDIMENSIONAL RADIAL POSITION, r/R 

Figure 29. Effect of Neglecting the Piper Cherokee 180 Fuselage Geometry 
(Replaced by a Short Afterbody) on the Computed Radial Distri- 
bution of Axial Velocity at the Propeller Plane 
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Figure 30a, Effect of Neglecting the Piper Cherokee 180 Wing on the 

Computed Radial Distribution of Tangential Velocity at the 
Propeller Plane 
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Figure 30b. Effect of Neglecting the Piper Cherokee 180 Wing on the' 
Computed Radial Distribution of Flow Angularity at the 
Propeller Plane (o * 2° Case) 
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Figtfre 30c. Effect of Neglecting the Piper Cherokee 180 Wing on the 
Computed Radial Distribution of Flow Angularity at the 
Propeller Plane (a * 10* Case) 
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Figure 31a, Effect of Wing Dihedral on the Computed Radial Distribution 
of Flow Angularity at the Propeller Plane of the Piper 
Cherokee 180 (o * 2° Case) 
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NONDIMEN SIONAL RADIAL POSITION, r/R 


Figure 31b. Effect of Wing Dihedral on the Computed Radial Distribution 
of Flow Angularity at the Propeller Plane of the Piper 
Cherokee 180 (a • 10° Case) 
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NONDIMENSIONAL RADIAL POSITION, r/R 


Figure 32a. Effect of Change In Cowl Inlet Inflow Velocity on the Computed 
Radial Distribution of Axial Velocity at the Propeller Plane 
of the Piper Cherokee 180 
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Figure 32b. Effect of Change in Cowl Inlet Inflow Velocity on the Computed 
Radial Distribution of Flow Angularity at the Propeller Platte 
of the Piper Cherokee 180 (a * 2° Case) 
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Figure 32c. Effect of Change in Cowl Inlet Inflow Velocity on the Computed 
Radial Distribution of Flow Angularity at the Propeller Plane 
of the Piper Cherokee 180 (a * 10° Case) 
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Figure 33a, Computed Azimuthal Distributions of Axial Velocity at the 
Propeller Plane of the Piper Cherokee 180 (at r/R ■ 0,35} 
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AZIMUTH POSITION, \ p (DEGREES) 


Figure 33h* Computed Azimuthal Distributions of Tangential Velocity- at the 
Propeller Plane of the Piper Cherokee 180 (at r/R - 0,35) 
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Figure 33c. Computed Azimuthal Distributions of Angle of Rotational Flow 
at the Propeller Plane of the Piper Cherokee 130 (at r/R ■ 0,35) 
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Figure 34a. Computed Azimuthal Distributions of Axial Velocity at the 
Propeller Plane of the Piper Cherokee 180 (at r/R * Q..5) 
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Figure 34b. Computed Azimuthal Distributions of Tangential Velocity at 
the Propeller Plane of the Piper Cherokee 180 (at r/R *.,0.,5) 
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Figure 35a. Computed Azimuthal Distributions of Axial Velocity- at the 
Propeller Plane of the Piper Cherokee 180 (at r/R ■ 0,75). 
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Figure 35b. Computed Azimuthal Distributions of Tangential Velocity- at 

the Propeller Plane of the Piper Cherokee 18Q Cat r/R « 0,75) 
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AZIMUTH POSITION, t|> (DEGREES) 

Figure 35c. Computed Azimuthal Distributions of Angle of Rotational Flow 
at the Propeller Plane of the Piper Cherokee 130 (at r/R " 0,75) 
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Figure 36a, Computed Azimuthal Distributions of Axial Velocity at the 
Propeller Plane of the Piper Cherokee 180 (at r/R - 1,0) 
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AZIMUTH POSITION, i|) (DEGREES) 

Figure 36b, Computed Azimuthal Distributions of Tangential Velocity at 
the Propeller Plane of the Piper Cherokee 180 Cat r/R ■ 1, 
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Figure 36c. Computed Azimuthal Distributions of Angle of Rotational flow 
at the Propeller Plane of the Piper Cherofcee 180 (at r/R«l,0) 
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Figure 37a. Computed Radial Distributions of Axial Velocity at the Pro- 
peller Plane of the Piper Cherokee 180 (at i|/ * 0®) 
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Figure 37b. Computed Radial Distributions of Tangential Velocity at the 
Propeller Plane of the Piper Cherokee 180 Cat ip ■ 0°) 
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NONDIMEN SIONAL RADIAL POSITION, r/R 

Figure 37c. Computed Radial Distributions of Angle of Rotational Flow at 
the Propeller Plane of the Piper Cherokee 18Q (at \p - 0°) 
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NONDIMENS IONAL RADIAL POSITION, r/R 

Figure 38a. Computed Radial Distributions of Axial Velocity at the 
Propeller Plane of the Piper Cherokee 180 Cat i p » 45 # ) 
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NONDIMENSIONAL RADIAL POSITION, r/R 

Figure 38h. Computed Radial Distributions of Tangential Velocity at the 
Propeller Plane of the Piper Cherokee 180 Cat m 45° 
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Figure 38c. Computed Radial Distributions of Angle of Rotational Flow at 
the Propeller Plane of the Piper Cherokee 180 Cat ip ■ 45 *) 
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NOND IMENSIONAL RADIAL POSITION, r/R 

Figure 39a, Computed Radial Distributions of Axial Velocity at the 
Propeller Plane of the Piper Cherokee 180 (at i|> - 90") 
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NONDIMENSIONAL RADIAL POSITION, r/R 


Figure 39b. Computed Radial Distributions of Tangential Velocity at the 
Propeller Plane of the Piper Cherokee 180 (at <p m 90°) 
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NOND IMEN S IONAL RADIAL POSITION, r/R 

Figure 39c. Computed Radial Distributions of Angle of Rotational Flow at 
the Propeller Plane of the Piper Cherokee 180 Cat ip ■ 90°) 
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NONDIMENSIONAL RADIAL POSITION, r/R 

Figure AOa. Computed Radial Distributions of Axial Velocity at the 
Propeller Plane of the Piper Cherokee 180 (at <Ji “ 135°) 
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NONDIMENSIONAL RADIAL POSITION, r/R 

Figure 40b. Computed Radial Distributions of Tangential Velocity at the 
Propeller Plane of the Piper Cherokee 180 (at ij> * 135°) 
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NONDIMEN SIONAL RADIAL POSITION, r/R 

. Computed Radial Distributions of Angle of Rotational Flow at 
the Propeller Plane of the Piper Cherokee 180 (at ■ 135°) 
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NONDIMENSIONAL RADIAL POSITION, r/R 

Figure 41a. Computed Radial Distributions of Axial Velocity at the 
Propeller Plane of the Piper Cherokee 180 (at ■ 180°) 
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NONDIMENS IONAL RADIAL POSITION, r/R 

Figure Alb. Computed Radial Distributions of Tangential Velocity at the 
Propeller Plane of the Piper Cherokee 180 (at iji “ 180°) 
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Figure 41c. Computed Radial Distributions of Angle of Rotational Flow at 
the Propeller Plane of the Piper Cherokee 180 (at ■ 180°) 
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NOND IMEN S IONAL RADIAL POSITION, r/R 

Figure 42a. Computed Radial Distributions of Axial Velocity at the 
Propeller Plane of the Piper Cherokee 180 (at ^ ” 225°) 
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Figure 42b. Computed Radial Distributions of Tangential Velocity at the 
Propeller Plane of the Piper Cherokee 180 Cat i|> - 225 e ) 
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Figure 42c. Computed Radial Distributions of Angle of Rotational Flow at 
the Propeller Plane of the Piper Cherokee 180 (at » 225°) 
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Computed Radial Distributions of Axial Velocity at the 
Propeller Plane of the Piper Cherokee 180 (at i ■ 270° 





TANGENTIAL VELOCITY RATIO 


ORIGINAL PAGE IS 
OF POOR QUALITY 


207 



NONDIMENSIONAL RADIAL POSITION, r/R 

Figure 43b. Computed Radial Distributions of Tangential Velocity at the 
Propeller Plane of the Piper Cherokee 180 Cat if) ■ 270°) 
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Figure 43c. Computed Radial Distributions of Angle of Rotational Flow at 
the Propeller Plane of the Piper Cherokee 180 (at ip " 270°) 
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NONDIMENSIONAL RADIAL POSITION, r/R 

Figure 44a. Computed Radial Distributions of Axial Velocity at the 
Propeller Plane of the Piper Cherokee 180 Cat ip • 315*) 
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NONDIMENS IONAL RADIAL POSITION, r/R 

Figure 44b. Computed Radial Distributions of Tangential Velocity at the 
Propeller Plane of the Piper Cherokee 180 (at - 315°) 
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NOND IMEN S IONAL RADIAL POSITION, r/R 

Figure 44c. Computed Radial Distributions of Angle of Rotational Flow at 
the Propeller Plane of the Piper Cherokee 180 (at ifi ■ 315°) 
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Figure A. la. Periphery Point Input on a Symmetric Body Cross Section 
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Figure A. lb. Periphery Point Input on a Nonsymmetric Body Cross Section 
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Figure A. 3. Quadrilateral Element With New Coplanar Corner Points 
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Figure A. 6a. Triangular Panel With Input Points 1 and 2 Coincident 
(Case a) 



Figure A. 6b. Triangular Panel With Input Points 3 and 4 Coincident 
(Case b) 
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Geometry and Induced Velocity of a Vortex Filament Having 
a Specified Locating and Orientation in Space 












